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Abstract 

We study Euler-Poincare systems (i.e., the Lagrangian analogue of Lie- 
Poisson Hamiltonian systems) defined on semidirect product Lie algebras. We 
first give a derivation of the Euler-Poincare equations for a parameter depen- 
dent Lagrangian by using a variational principle of Lagrange d'Alembert type. 
Then we derive an abstract Kelvin-Noether theorem for these equations. We 
also explore their relation with the theory of Lie-Poisson Hamiltonian systems 
defined on the dual of a semidirect product Lie algebra. The Legendre trans- 
formation in such cases is often not invertible; so it does not produce a corre- 
sponding Euler-Poincare system on that Lie algebra. We avoid this potential 
difficulty by developing the theory of Euler-Poincare systems entirely within 
the Lagrangian framework. We apply the general theory to a number of known 
examples, including the heavy top, ideal compressible fluids and MHD. We 
also use this framework to derive higher dimensional Camassa-Holm equations, 
which have many potentially interesting analytical properties. These equations 
are Euler-Poincare equations for geodesies on diffeomorphism groups (in the 
sense of the Arnold program) but where the metric is H 1 rather than L 2 . 
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1 History and Background 

To put our paper in context, we shall pick up the thread of the history of mechanics 
in the later part of the 1800's. By that time, through the work of many people, 
including Euler, Lagrange, Hamilton, Jacobi, and Routh, it was well understood 
that the equations of mechanics are expressible in either Hamiltonian or Lagrangian 
form. 

The Lagrangian formulation of mechanics can be based on the variational prin- 
ciples behind Newton's fundamental laws of force balance F = ma. One chooses a 
configuration space Q (a manifold, assumed to be of finite dimension n to start the 
discussion) with coordinates denoted q l ,i = 1,... , n, that describe the configura- 
tion of the system under study. One then forms the velocity phase space TQ (the 
tangent bundle of Q). Coordinates on TQ are denoted (q , . . . , q n , q , . . . , q n ), and 
the Lagrangian is regarded as a function L : TQ — ► R. In coordinates, one writes 
L(q l , q l ,t), which is shorthand notation for L(q , . . . , q n , q , q n ,t). Usually, L is 
the kinetic minus the potential energy of the system and one takes q % = dq l /dt to be 
the system velocity. The variational principle of Hamilton states that the variation 
of the action is stationary at a solution: 

5& = 5 I L{q\q\t)dt = 0. (1.1) 

J a 

In this principle, one chooses curves q l {t) joining two fixed points in Q over a fixed 
time interval [a, b], and calculates the action 6, which is the time integral of the 
Lagrangian, regarded as a function of this curve. Hamilton's principle states that 
the action 6 has a critical point at a solution in the space of curves. As is well 



2 



known, Hamilton's principle is equivalent to the Euler-Lagrange equations: 

d dL dL . 

If the system is subjected to external forces, these are to be added to the right 
hand side of the Euler-Lagrange equations. For the case in which L comprises kinetic 
minus potential energy, the Euler-Lagrange equations reduce to a geometric form of 
Newton's second law. For Lagrangians that are purely kinetic energy, it was already 
known in Poincare's time that the corresponding solutions of the Euler-Lagrange 
equations are geodesies. (This fact was certainly known to Jacobi by 1840, for 
example.) 

To pass to the Hamiltonian formalism, one introduces the conjugate momenta 

dL 

Pi = W i = 1 '---' n ' 

and makes the change of variables (g\(f) i— > (q l ,Pi), by a Legendre transformation. 
The Lagrangian is called regular when this change of variables is invertible. The 
Legendre transformation introduces the Hamiltonian 

n 

H(q\pi,t) = J2PjQ 3 ~ L{q\q\t). (1.4) 

One shows that the Euler-Lagrange equations are equivalent to Hamilton's equa- 
tions: 

dq^ = 3H_ dm _ (M 

dt dpi 1 dt dq v [ ' ' 

where i = 1, . . . ,n. There are analogous Hamiltonian partial differential equations 
for field theories such as Maxwell's equations and the equations of fluid and solid 
mechanics. 

Hamilton's equations can be recast in Poisson bracket form as 

F = {F,H}, (1.6) 
where the canonical Poisson brackets are given by 

{F,G} = ^ (—— -— —A . (1.7) 

Associated to any configuration space Q is a phase space T*Q called the cotan- 
gent bundle of Q, which has coordinates (q 1 , . . . , q n ,pi, ■ ■ ■ ,p n )- On this space, the 
canonical Poisson bracket is intrinsically defined in the sense that the value of {F, G} 
is independent of the choice of coordinates. Because the Poisson bracket satisfies 
{F, G} = —{G, F} and in particular {H, H} = 0, we see that H = 0; that is, energy 
is conserved along solutions of Hamilton's equations. This is the most elementary 
of many deep and beautiful conservation properties of mechanical systems. 
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Poincare and the Euler equations. Poincare played an enormous role in the 
topics treated in the present paper. We mention a few of Poincare's contributions 
that are relevant here. First is his work on the gravitating fluid problem, continuing 
the line of investigation begun by MacLaurin, Jacobi and Riemann. Some solutions 
of this problem still bear his name today. This work is summarized in Chandrasekhar 
[1967, 1977] (see Poincare [1885, 1890, 1892, 1901a] for the original treatments). This 
background led to his famous paper, Poincare [1901b], in which he laid out the basic 
equations of Euler type, including the rigid body, heavy top and fluids as special 
cases. Abstractly, these equations are determined once one is given a Lagrangian on 
a Lie algebra. We shall make some additional historical comments on this situation 
below, after we present a few more mechanical preliminaries. It is because of the 
paper Poincare [1901b] that the name Euler Poincare equations is now used for 
these equations. 

To state the Euler-Poincare equations, let q be a given Lie algebra and let 
/ : g — > R be a given function (a Lagrangian), let £ be a point in q and let / 6 g* 
be given forces (whose nature we shall explicate later). Then the evolution of the 
variable £ is determined by the Euler-Poincare equations. Namely, 



The notation is as follows: dl/d^ G 0* (the dual vector space) is the derivative of I 
with respect to £; we use partial derivative notation because I is a function of the 
vector £ and because shortly I will be a function of other variables as well. The map 
ad^ : q — ► q is the linear map r] i-> [£,tj], where [£,77] denotes the Lie bracket of £ 
and 77, and where ad| : 5* — > q* is its dual (transpose) as a linear map. In the case 
that / = 0, we will call these equations the basic Euler Poincare equations. 

These equations are valid for either finite or infinite dimensional Lie algebras. For 
fluids, Poincare was aware that one needs to use infinite dimensional Lie algebras, 
as is clear in his paper Poincare [1910]. He was aware that one has to be careful 
with the signs in the equations; for example, for rigid body dynamics one uses the 
equations as they stand, but for fluids, one needs to be careful about the conventions 
for the Lie algebra operation ad^; cf. Chetayev [1941]. 

To state the equations in the finite dimensional case in coordinates, one must 
choose a basis e\, . . . , e r of q (so dimg = r). Define, as usual, the structure constants 
C^ b of the Lie algebra by 



where a, b run from 1 to r. If £ G 0, its components relative to this basis are 
denoted £ a . If e 1 , . . . , e n is the corresponding dual basis, then the components of 
the differential of the Lagrangian I are the partial derivatives dl/d£ a . The Euler- 
Poincare equations in this basis are 



r 




(1.8) 



d=i 




(1.9) 
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For example, consider the Lie algebra M 3 with the usual vector cross product. (Of 
course, this is the Lie algebra of the proper rotation group in M 3 .) For I : M 3 — ► R, 
the Euler-Poincare equations become 

dt on on x 

which generalize the Euler equations for rigid body motion. 

These equations were written down for a certain class of Lagrangians I by La- 
grange [1788, Volume 2, Equation A on p. 212], while it was Poincare [1901b] who 
generalized them (without reference to the ungeometric Lagrange!) to an arbitrary 
Lie algebra. However, it was Lagrange who was grappeling with the derivation and 
deeper understanding of the nature of these equations. While Poincare may have 
understood how to derive them from other principles, he did not reveal this. 

Of course, there was a lot of mechanics going on in the decades leading up to 
Poincare 's work and we shall comment on some of it below. However, it is a curious 
historical fact that the Euler-Poincare equations were not pursued extensively until 
quite recently. While many authors mentioned these equations and even tried to 
understand them more deeply (see, e.g., Hamel [1904, 1949] and Chetayev [1941]), it 
was not until the Arnold school that this understanding was at least partly achieved 
(see Arnold [1966a, c] and Arnold [1988]) and was used for diagnosing hydrodynam- 
ical stability (e.g., Arnold [1966b]). 

It was already clear in the last century that certain mechanical systems resist 
the usual canonical formalism, either Hamiltonian or Lagrangian, outlined in the 
first paragraph. The rigid body provides an elementary example of this. In another 
example, to obtain a Hamiltonian description for ideal fluids, Clebsch [1857, 1859] 
found it necessary to introduce certain nonphysical potentials^. 



More about the rigid body. In the 

equations are usually written as follows: 

m 2 = (h 

h^3 = {h 



absence of external forces, the rigid body 
-Wi, (1-10) 

-I 2 )fil^2, 



where fi = (Oi,^,^) is the body angular velocity vector and I\,l2,l3 are the 
moments of inertia of the rigid body. Are these equations as written Lagrangian or 
Hamiltonian in any sense? Since there are an odd number of equations, they cannot 
be put in canonical Hamiltonian form. 

One answer is to reformulate the equations on TSO(3) or T*SO(3), as is clas- 
sically done in terms of Euler angles and their velocities or conjugate momenta, 
relative to which the equations are in Euler-Lagrange or canonical Hamiltonian 

1 For modern accounts of Clebsch potentials and further references, see Holm and Kupersh- 
midt [1983], Marsden and Weinstein [1983], Marsden, Ratiu, and Weinstein [1984a, b], Cendra and 
Marsden [1987], Cendra, Ibort, and Marsden [1987] and Goncharov and Pavlov [1997]. 
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form. However, this reformulation answers a different question for a six dimensional 
system. We are interested in these structures for the equations as given above. 

The Lagrangian answer is easy: these equations have Euler-Poincare form on 
the Lie algebra M 3 using the Lagrangian 



j(n) = i(iifi? + j 2 fi! + j 3 fi§). 



(l.ii) 



which is the (rotational) kinetic energy of the rigid body. 

One of our main messages is that the Euler-Poincare equations possess a natural 
variational principle. In fact, the Euler rigid body equations are equivalent to the 
rigid body action principle 



56 red = 5 ldt = 0, 

J a 

where variations of fi are restricted to be of the form 

^n = s + nxs, 



(1.12) 



(1.13) 



in which S is a curve in M 3 that vanishes at the endpoints. As before, we regard 
the reduced action S re( j as a function on the space of curves, but only consider 
variations of the form described. The equivalence of the rigid body equations and 
the rigid body action principle may be proved in the same way as one proves that 
Hamilton's principle is equivalent to the Euler-Lagrange equations: Since I (ft) = 
|(Ifi,f2), and I is symmetric, we obtain 



5 Idt 



(ISl,SSl)dt 



(in, s + si x s) dt 



in,E) + (in,n x s) 



d. 
~dl l 



b 1 d 

—Ifl + m x ft, £ ) dt, 

dt 



upon integrating by parts and using the endpoint conditions, £(&) 
Since 5] is otherwise arbitrary, (|1.12 ) is equivalent to 



Sfa) 



-— (m) +in x ft = o, 

which are Euler's equations. 

Let us explain in concrete terms (that will be abstracted later) how to derive 
this variational principle from the standard variational principle of Hamilton. 

We regard an element R € SO(3) giving the configuration of the body as a map 
of a reference configuration B C I 3 to the current configuration R(S); the map R 
takes a reference or label point X S B to a current point x = R(^") £ R(S). When 



6 



the rigid body is in motion, the matrix R is time-dependent and the velocity of a 
point of the body is x = RX = RR _1 x. Since R is an orthogonal matrix, R _1 R 
and RR 1 are skew matrices, and so we can write 

x = RR _1 x = u x x, (1-14) 

which defines the spatial angular velocity vector u. Thus, u is essentially given 
by right translation of R to the identity. 

The corresponding body angular velocity is defined by 

fi = R "V (1.15) 

so that f2 is the angular velocity relative to a body fixed frame. Notice that 

R ^RI = R X RR ~ l x = R-^u; x x) 

= R ^ x BT l x = fixX, (1.16) 

so that fi is given by left translation of R to the identity. The kinetic energy is 
obtained by summing up m\x\ 2 /2 (where |»| denotes the Euclidean norm) over the 
body: 

K = \ [ p(X)\RX\ 2 d 3 X, (1.17) 



2JB 

in which p is a given mass density in the reference configuration. Since 

|RX| = \u>xx\ = |R _1 (^ x x)\ = \flx X\, 
K is a quadratic function of CI. Writing 

K = htflQ, (1.18) 

defines the moment of inertia tensor I, which, provided the body does not 
degenerate to a line, is a positive-definite (3 x 3) matrix, or better, a quadratic form. 
This quadratic form can be diagonalized by a change of basis; thereby defining the 
principal axes and moments of inertia. In this basis, we write I = diag(ii, 12,1s). The 
function K is taken to be the Lagrangian of the system on TSO(3) (and by means of 
the Legendre transformation we obtain the corresponding Hamiltonian description 



on T*SO(3)). Notice that K in equation (1.17) is left (not right) invariant on 
TSO(3). It follows that the corresponding Hamiltonian is also left invariant. 

In the Lagrangian framework, the relation between motion in R space and mo- 
tion in body angular velocity (or ft) space is as follows: The curve R(t) G SO(3) 
satisfies the Euler-Lagrange equations for 

L(R,R) = - [ p(X)\iiX\ 2 d 3 X, (1.19) 
2 Jb 
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if and only if Q(t) denned by R 1 Rv = fi x v for all v E IR 3 satisfies Euler's 
equations 

in = lO X ft. (1.20) 

An instructive proof of this relation involves understanding how to reduce vari- 
ational principles using their symmetry groups. By Hamilton's principle, R(i) sat- 
isfies the Euler-Lagrange equations, if and only if 

5 I L dt = 0. 



Let l(Q) = I (in) • n, so that I (ft) = L(R, R) if R and n are related as above. 
To see how we should transform Hamilton's principle, define the skew matrix n by 
Civ = n x v for any v E M 3 , and differentiate the relation R _1 R = Cl with respect 
to R to get 

-R-^^R^R + R^OSR) = 6TI. (1.21) 
Let the skew matrix E be defined by 

E = R _1 <5R, (1.22) 

and define the vector E by 

Sv = E x v. (1.23) 

Note that 

S = -R X RR X 5R + R 1 JR, 

so 

R _1 5R = E + R^RS . (1.24) 
Substituting flL2§ and ( pT22| ) into (|L2l] ) gives 

-En + E + nE = sn, 

that is, 

jn = E + [n,E]. (i.25) 

The identity [n, S] = (O X E)" holds by Jacobi's identity for the cross product and 
so 

5n = E + nxE. (1.26) 
These calculations prove the following: 



8 



Theorem 1.1 Hamilton's variational principle 

56 = 5 f Ldt = (1.27) 

J a 

on TSO(3) is equivalent to the reduced variational principle 

5& Tcd = 5 f ldt = (1.28) 

J a 

on M 3 where the variations 5CI are of the form (1.26) with E(a) = E(6) = 0. 



Hamiltonian Form. If, instead of variational principles, we concentrate on Pois- 
son brackets and drop the requirement that they be in the canonical form, then there 
is also a simple and beautiful Hamiltonian structure for the rigid body equations 
that is now well knownQ. To recall this, introduce the angular momenta 

dL 

11; = /;^ = —, 1 = 1,2,3, (1.29) 
so that the Euler equations become 

rii = ^r^n 2 n 3 , 

ri 2 = ^— — n^, (i.30) 

that is, 

ri = nxfi. (1.31) 

Introduce the following rigid body Poisson bracket on functions of the ITs: 

{F, G}(n) = -n • (VnF x V n G) (1.32) 

and the Hamiltonian 

One checks that Euler's equations are equivalent to F = {F,H}. 

The rigid body variational principle and the rigid body Poisson bracket are 
special cases of general constructions associated to any Lie algebra q. Since we have 
already described the general Euler-Poincare construction on q, we turn next to the 
Hamiltonian counterpart on the dual space. 



2 See Marsden and Ratiu [1994] for details, references, and the history of this structure. 
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The Lie-Poisson Equations. Let F, G be real valued functions on the dual space 
g*. Denoting elements of g* by /j,, let the functional derivative of F at /x be the unique 
element 5F/Sfi of g defined by 

lim ^[Fiji + e5^) - Fiji)] = (6^, ^ , (1.34) 

for all Sfi € g*, where (,) denotes the pairing between g* and g. Define the (±) 
Lie-Poisson brackets by 



{F,G} ± ( M )=±(/i 



£G' 



(1.35) 



Using the coordinate notation introduced above, the (±) Lie-Poisson brackets be- 
come 

where // = YJd=i VdZ d ■ 

The Lie-Poisson equations, determined by F = {F, H} read 

L v^ ^ OH 

b,d=l 0fib 

or intrinsically, 

V = T&d* dH/dfif i. (1.37) 

This setting of mechanics is a special case of the general theory of systems on 
Poisson manifolds, for which there is now an extensive theoretical development. 
(See Guillemin and Sternberg [1984] and Marsden and Ratiu [1994] for a start on 
this literature.) There is an especially important feature of the rigid body bracket 
that carries over to general Lie algebras, namely, Lie-Poisson brackets arise from 
canonical brackets on the cotangent bundle (phase space) T*G associated with a Lie 
group G which has g as its associated Lie algebra. 

For a rigid body which is free to rotate about its center of mass, G is the 
(proper) rotation group SO (3). The choice of T*G as the primitive phase space is 
made according to the classical procedures of mechanics described earlier. For the 
description using Lagrangian mechanics, one forms the velocity-phase space TSO(3). 
The Hamiltonian description on T*G is then obtained by standard procedures. 

The passage from T*G to the space of II's (body angular momentum space) 
is determined by left translation on the group. This mapping is an example of a 
momentum map; that is, a mapping whose components are the "Noether quantities" 
associated with a symmetry group. The map from T*G to g* being a Poisson 
(canonical) map is a general fact about momentum maps. The Hamiltonian point 
of view of all this is again a well developed subject. 
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Geodesic motion. As emphasized by Arnold [1966a], in many interesting cases, 
the Euler-Poincare equations on a Lie algebra g correspond to geodesic motion on 
the corresponding group G. We shall explain the relationship between the equations 



on q and on G shortly, in theorem [L2 . Similarly, on the Hamiltonian side, the 
preceding paragraphs explained the relation between the Hamiltonian equations on 
T*G and the Lie-Poisson equations on q*. However, the issue of geodesic motion is 
simple: if the Lagrangian or Hamiltonian on q or q* is purely quadratic, then the 
corresponding motion on the group is geodesic motion. 



More History. The Lie-Poisson bracket was discovered by Sophus Lie (Lie [1890], 
Vol. II, p. 237). However, Lie's bracket and his related work was not given much 
attention until the work of Kirillov, Kostant, and Souriau (and others) revived it in 
the mid-1960s. Meanwhile, it was noticed by Pauli and Martin around 1950 that 
the rigid body equations are in Hamiltonian form using the rigid body bracket, but 
they were apparently unaware of the underlying Lie theory. It would seem that 
while Poincare was aware of Lie theory, in his work on the Euler equations he was 
unaware of Lie's work on Lie-Poisson structures. He also seems not to have been be 
aware of the variational structure of the Euler equations. 



The heavy top. Another system important to Poincare and also for us in this 
paper is the heavy top; that is, a rigid body with a fixed point in a gravitational field. 
For the Lie-Poisson description, the underlying Lie algebra, surprisingly, consists of 
the algebra of infinitesimal Euclidean motions in M 3 . These do not arise as actual 
Euclidean motions of the body since the body has a fixed point! As we shall see, 
there is a close parallel with the Poisson structure for compressible fluids. 

The basic phase space we start with is again T*SO(3). In this space, the equa- 
tions are in canonical Hamiltonian form. Gravity breaks the symmetry and the 
system is no longer SO(3) invariant, so it cannot be written entirely in terms of the 
body angular momentum II. One also needs to keep track of T, the "direction of 
gravity" as seen from the body (T = R _1 k where the unit vector k points upward 
and R is the element of SO(3) describing the current configuration of the body). 
The equations of motion are 

rii = ^-^n 2 n 3 + Mgi (rV - rV), 

ri 2 = ^— ^n 3 n! + Mgi (rV - rV), (i.38) 

ri 3 = ^f^n^ + Mgi (rV - rV), 



or, in vector notation, 



and 



tl = Uxn + Mg£Txx, (1-39) 

r = r x n, (i.4o) 



n 



where M is the body's mass, g is the acceleration of gravity, x is t ne unit vector 
on the line connecting the fixed point with the body's center of mass, and I is the 
length of this segment. 

The Lie algebra of the Euclidean group is se(3) = R 3 x I 3 with the Lie bracket 

[(£, u), (r?, v)] = {i x rj, £ x v - rj x u). (1.41) 

We identify the dual space with pairs (II, T); the corresponding (— ) Lie-Poisson 
bracket called the heavy top bracket is 

{F,G}(n,T) = -n-(V n Fx V n G) 

-r.(V n FxVrG-VnGxV r F). (1.42) 

The above equations for n, F can be checked to be equivalent to 

F = {F,H}, (1.43) 

where the heavy top Hamiltonian 

i /tt2 tt 2 TT 2 \ 

H(U, r) = -(-^ + ^ + ^)+Mg£T. X (1-44) 

is the total energy of the body (see, for example, Sudarshan and Mukunda [1974]). 

The Lie algebra of the Euclidean group has a structure which is a special case of 
what is called a semidirect product. Here it is the product of the group of rotations 
with the translation group. It turns out that semidirect products occur under rather 
general circumstances when the symmetry in T*G is broken. In particular, there 
are similarities in structure between the Poisson bracket for compressible flow and 
that for the heavy top. The general theory for semidirect products will be reviewed 
shortly. 

A Kaluza-Klein form for the heavy top. We make a remark about the heavy 
top equations that is relevant for later purposes. Namely, since the equations have a 
Hamiltonian that is of the form kinetic plus potential, it is clear that the equations 
are not of Lie-Poisson form on so(3)* ; the dual of the Lie algebra of SO(3) and 
correspondingly, are not geodesic equations on SO (3). While the equations are Lie- 
Poisson on se(3)*, the Hamiltonian is not quadratic, so again the equations are not 
geodesic equations on SE(3). 

However, they can be viewed a different way so that they become Lie-Poisson 
equations for a different group and with a quadratic Hamiltonian. In particular, they 
are the reduction of geodesic motion. To effect this, one changes the Lie algebra from 
se(3) to the product se(3) x so (3). The dual variables are now denoted II, T, x- We 
regard the variable % as a momentum conjugate to a new variable, namely a ghost 
element of the rotation group in such a way that X is a constant of the motion; in 
Kaluza-Klein theory for charged particles on thinks of the charge this way, as being 
the momentum conjugate to a (ghost) cyclic variable. 
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We modify the Hamiltonian by replacing T-x by f° r example, r-x+||r|| 2 + ||x|| 2 , 
or any other terms of this sort that convert the potential energy into a positive 
definite quadratic form in T and X- The added terms, being Casimir functions, 
do not affect the equations of motion. However, now the Hamiltonian is purely 
quadratic and hence comes from geodesic motion on the group SE(3) x SO(3). 
Notice that this construction is quite different from that of the well known Jacobi 
metric method. 

Later on in our study of continuum mechanics, we shall repeat this construction 
to achieve geodesic form for some other interesting continuum models. Of course 
one can also treat a heavy top that is charged or has a magnetic moment using these 
ideas. 



Incompressible Fluids. Arnold [1966a] showed that the Euler equations for an 
incompressible fluid could be given a Lagrangian and Hamiltonian description simi- 
lar to that for the rigid body. His approach^ has the appealing feature that one sets 
things up just the way Lagrange and Hamilton would have done: one begins with 
a configuration space Q, forms a Lagrangian L on the velocity phase space TQ and 
then Legendre transforms to a Hamiltonian H on the momentum phase space T*Q. 
Thus, one automatically has variational principles, etc. For ideal fluids, Q = G is 
the group Diff vo i(2?) of volume preserving transformations of the fluid container (a 
region D in IR 2 or IR 3 , or a Riemannian manifold in general, possibly with boundary). 
Group multiplication in G is composition. 

The reason we select G = Diff vo i(I?) as the configuration space is similar to that 
for the rigid body; namely, each tp in G is a mapping of T> to T> which takes a 
reference point X G V to a current point x = f{X) G V; thus, knowing <p tells us 
where each particle of fluid goes and hence gives us the current fluid configuration . 
We ask that i^bea diffeomorphism to exclude discontinuities, cavitation, and fluid 
interpenetration, and we ask that tp be volume preserving to correspond to the 
assumption of incompressibility. 

A motion of a fluid is a family of time-dependent elements of G, which we write 
as x = ip(X,t). The material velocity field is defined by V(X, t) = dcp(X,t)/dt, 
and the spatial velocity field is defined by v(x,t) = ~V(X,t) where x and X are 
related by x = (p(X,t). If we suppress "i" and write tp for V, note that 

v = l po(p- 1 i.e., V( = V[0^" 1 , (1-45) 



where <pt(x) = (p(X,t). We can regard ( 1.45| ) as a map from the space of (<p,<p) 



(material or Lagrangian description) to the space of v's (spatial or Eulerian descrip- 
tion). Like the rigid body, the material to spatial map ( |1.45| ) takes the canonical 
bracket to a Lie-Poisson bracket; one of our goals is to understand this reduction. 
Notice that if we replace <p> by <p o r/ for a fixed (time- independent) r/ G Diff vo i(D), 
then Lp o ip^ 1 is independent of rj; this reflects the right invariance of the Eulerian 

3 Arnold's approach is consistent with what appears in the thesis of Ehrenfest from around 
1904; see Klein [1970]. However, Ehrenfest bases his principles on the more sophisticated curvature 
principles of Gauss and Hertz. 
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description (v is invariant under composition of ip by r/ on the right). This is also 
called the particle relabeling symmetry of fluid dynamics. The spaces TG and 
T*G represent the Lagrangian (material) description and we pass to the Eulerian 
(spatial) description by right translations and use the (+) Lie-Poisson bracket. One 
of the things we shall explain later is the reason for the switch between right and 
left in going from the rigid body to fluids. 

The Euler equations for an ideal, incompressible, homogeneous fluid moving 
in the region T> are 

_ + (v • V)v = -Vp (1.46) 

with the constraint div v = and boundary conditions: v is tangent to &D. 

The pressure p is determined implicitly by the divergence-free (volume preserv- 
ing) constraint div v = 0. The associated Lie algebra q is the space of all divergence- 
free vector fields tangent to the boundary. This Lie algebra is endowed with the 
negative Jacobi-Lie bracket of vector fields given by 

3=1 

(The subscript L on [• , •] refers to the fact that it is the left Lie algebra bracket on 
g. The most common convention for the Jacobi-Lie bracket of vector fields, also the 
one we adopt, has the opposite sign.) We identify g and q* by using the pairing 

(v,w) = / v wd 3 x. (1.48) 
Jv 

Hamiltonian structure for fluids. Introduce the (+) Lie-Poisson bracket, called 
the ideal fluid bracket, on functions of v by 

{F,G}(v)= / v \^f,f] d 3 x, (1.49) 

where 5F/5v is defined by 

lim-[F(v + etfv)-F(v)l = / f 5v • ^ ] d 3 x. (1.50) 

With the energy function chosen to be the kinetic energy, 

H(v) = l I |v| 2 d 3 x, (1.51) 



one can verify that the Euler equations (1.46) are equivalent to the Poisson bracket 
equations 

F = {F,H} (1.52) 
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for all functions F on q. For this, one uses the orthogonal decomposition w = 
Pw + Vp of a vector field w into a divergence- free part Pw in q and a gradient. The 
Euler equations can be written as 



One can also express the Hamiltonian structure in terms of the vorticity as a 
basic dynamic variable and show that the preservation of coadjoint orbits amounts to 
Kelvin's circulation theorem. We shall see a Lagrangian version of this property later 
in the paper. Marsden and Weinstein [1983] show that the Hamiltonian structure 
in terms of Clebsch potentials fits naturally into this Lie-Poisson scheme, and that 
Kirchhoff 's Hamiltonian description of point vortex dynamics, vortex filaments, and 
vortex patches can be derived in a natural way from the Hamiltonian structure 
described above. 

Lagrangian structure for fluids. The general framework of the Euler-Poincare 
and the Lie-Poisson equations gives other insights as well. For example, this general 
theory shows that the Euler equations are derivable from the "variational principle" 



where u is a vector field (representing the infinitesimal particle displacement) van- 
ishing at the temporal endpoints. The constraints on the allowed variations of the 
fluid velocity field are commonly known as "Lin constraints" and their nature was 
clarified by Newcomb [1962] and Bretherton [1970]. This itself has an interesting 
history, going back to Ehrenfest, Boltzmann, and Clebsch, but again, there was little 
if any contact with the heritage of Lie and Poincare on the subject. 

The Basic Euler— Poincare Equations. We now recall the abstract derivation 
of the "basic" Euler-Poincare equations (i.e., the Euler-Poincare equations with no 
forcing or advected parameters) for left-invariant Lagrangians on Lie groups (see 
Marsden and Scheurle [1993a,b], Marsden and Ratiu [1994] and Bloch et al. [1996]). 

Theorem 1.2 Let G be a Lie group and L : TG — ► M a left (respectively, right) 
invariant Lagrangian. Let I : q — > R be its restriction to the tangent space at the 
identity. For a curve g(t) G G, let = g(t)- l g(t); i.e., = T g ^L g ^-ig(t) 
(respectively, £(t) = g{t)g{t)~ l ). Then the following are equivalent: 

i Hamilton's principle 



— + P(v ■ Vv) = 0. 



(1.53) 




which should hold for all variations 5v of the form 



<5v = u + [u, v]l 




(1.54) 



holds, as usual, for variations Sg(t) of g(t) vanishing at the endpoints. 
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ii The curve g{t) satisfies the Euler- Lagrange equations for L on G. 
iii The "variational" principle 

5 I l(£(t))dt = (1.55) 

J a 

holds on q, using variations of the form 

5Z = ri±[ti,v], (1-56) 

where n vanishes at the endpoints (+ corresponds to left invariance and — to 
right invariance) F\ 



iv The basic Euler Poincare equations hold 

dSl 51 

Jt5i- ±ad zTz 



±ad| — . (1.57) 



Basic Ideas of the Proof. First of all, the equivalence of i and ii holds on the 
tangent bundle of any configuration manifold Q, by the general Hamilton principle. 
To see that ii and iv are equivalent, one needs to compute the variations <5£ induced 
on £ = g^g = TL g ~ig by a variation of g. We will do this for matrix groups; see 
Bloch, Krishnaprasad, Marsden, and Ratiu [1994] for the general case. To calculate 
this, we need to differentiate g~ l g in the direction of a variation 5g. If 5g = dg/de 
at e = 0, where g is extended to a curve g e , then, 



while if rj = g 1 5g, then 



kc d -i d 
5i = Je 9 dt 9 > 



. _ d _2 d 
^ dt 9 de 9 



The difference <5£ — r) is thus the commutator [£,??]. 

To complete the proof, we show the equivalence of iii and iv in the left-invariant 
case. Indeed, using the definitions and integrating by parts produces, 



5 J l(0dt = j ^ dt = f Jz(rj + ^rj)dt 

n dt , 



d ( 5l\ ^51 



dt\5a +ad( si 

so the result follows. 

There is of course a right invariant version of this theorem in which £ = gg~ 1 
and the Euler-Poincare equations acquire appropriate minus signs as in equation 
( 1.57j ). We shall go into this in detail later. 

4 Because there are constraints on the variations, this principle is more like a Lagrange d'Alem- 
bert principle, which is why we put "variational" in quotes. As we shall explain, such problems are 
not literally variational. 
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Since the Euler-Lagrange and Hamilton equations on TQ and T*Q are equivalent 
in the regular case, it follows that the Lie-Poisson and Euler-Poincare equations 
are then also equivalent. To see this directly, we make the following Legendre 
transformation from q to g*: 



Note that 



Sh I bi\ /SI 6£\ ( 



and so it is now clear that the Lie-Poisson equations ( 1.37| ) and the Euler-Poincare 
equations (|l.57| ) are equivalent. 

We close this paragraph by mentioning the geodesic property of the basic Euler- 
Poincare form. When I is a metric on TG, the basic Euler-Poincare equations are the 
geodesic spray equations for geodesic motion on the group G with respect to that 
metric. For discussions of this property in applications, see, e.g., Arnold [1966a] 
for the Euler equations of an incompressible ideal fluid, and Ovsienko and Khesin 
[1987] for the KdV shallow water equation. (An account of the latter case from the 
Euler-Poincare viewpoint may also be found in Marsden and Ratiu [1994].) Zeitlin 
and Pasmanter [1994] discuss the geodesic property for certain ideal geophysical 
fluid flows; Zeitlin and Kambe [1993] and Ono [1995a, 1995b] discuss it for ideal 
MHD; and Kouranbaeva [1997] for the integrable Camassa-Holm equation. From 
one viewpoint, casting these systems into basic Euler-Poincare form explains why 
they share the geodesic property. 



Lie-Poisson Systems on Semidirect Products. As we described above, the 
heavy top is a basic example of a Lie-Poisson Hamiltonian system defined on the 
dual of a semidirect product Lie algebra. The general study of Lie-Poisson equations 
for systems on the dual of a semidirect product Lie algebra grew out of the work 
of many authors including Sudarshan and Mukunda [1974], Vinogradov and Ku- 
pershmidt [1977], Ratiu [1980], Guillemin and Sternberg [1980], Ratiu [1981, 1982], 
Marsden [1982], Marsden, Weinstein, Ratiu, Schmidt and Spencer [1983], Holm and 
Kupershmidt [1983], Kupershmidt and Ratiu [1983], Holmes and Marsden [1983], 
Marsden, Ratiu and Weinstein [1984a, b], Guillemin and Sternberg [1984], Holm, 
Marsden, Ratiu and Weinstein [1985], Abarbanel, Holm, Marsden, and Ratiu [1986] 
and Marsden, Misiolek, Perlmutter and Ratiu [1997]. As these and related refer- 
ences show, the Lie-Poisson equations apply to a wide variety of systems such as 
the heavy top, compressible flow, stratified incompressible flow, and MHD (magne- 
tohydrodynamics) . We review this theory in §2 below. 

In each of the above examples as well as in the general theory, one can view the 
given Hamiltonian in the material representation as one that depends on a param- 
eter; this parameter becomes dynamic when reduction is performed; this reduction 
amounts in many examples to expressing the system in the spatial representation. 
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Goals of this Paper. The first goal of this paper is to study a Lagrangian ana- 
logue of the Hamiltonian semidirect product theory. The idea is to carry out a 
reduction for a Lagrangian that depends on a parameter and to use the ideas of 
reduction of variational principles from Marsden and Scheurle [1993a, b] and Bloch, 
Krishnaprasad, Marsden and Ratiu [1996] to directly reduce the problem to one 
that parallels Lie-Poisson systems on the duals of semidirect products. We call 
the resulting equations the Euler-Poincare equations since, as we have explained, 
Poincare [1901b] came rather close to this general picture. These equations gener- 
alize the basic Euler-Poincare equations on a Lie algebra in that they depend on a 
parameter and this parameter in examples has the interpretation of being advected, 
or Lie dragged, as is the density in compressible flow. 

One of the reasons this process is interesting and cannot be derived directly 
from its Hamiltonian counterpart by means of the Legendre transformation is that 
in many examples, such as the heavy top, the Hamiltonian describing the Lie-Poisson 
dynamics is degenerate; that is, the Legendre transformation is not invertible. 

A second major goal is to prove a version of the Noether theorem in an action 
principle formulation that leads immediately to a Kelvin circulation type theorem for 
continuum mechanics. We call this general formulation the Kelvin-Noether theorem. 

Finally, we provide a number of applications of the Euler-Poincare equations in 
ideal continuum dynamics which illustrate the power of this approach in unifying 
various known models, as well as in formulating new models. We also discuss some 
circumstances when the equations can be cast into the form of geodesies on certain 
infinite dimensional groups. 



Outline of the remainder of this paper. In the next section we review the 
semidirect product theory for Hamiltonian systems. Then in section || we consider 
the Lagrangian counterpart to this theory. Section 4] discusses the Kelvin-Noether 
theorem for the Euler-Poincare equations. Section | illustrates the general theory 
in the example of the heavy top. We introduce the Euler-Poincare equations for 
continua in section |6| and consider their applications to compressible flow (including 
MHD and adiabatic Maxwell-fluid plasmas) in section [|. Various approximate forms 
of the shallow water equations, such as the Boussinesq equations, the Camassa-Holm 
equation and its new higher-dimensional variants are developed in section ||. In other 
publications, the Maxwell- Vlasov equations will be considered as well as a general 
framework for the theory of reduction by stages. 

In the remainder of this paper we assume that the reader is familiar with Lie- 
Poisson Hamiltonian systems defined on duals of Lie algebras and the Lie-Poisson 
reduction theorem, reviewed above. We refer to Marsden and Ratiu [1994] for a 
detailed exposition of these matters. 
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2 Hamiltonian Semidirect Product Theory 



We first recall how the Hamiltonian theory proceeds for systems defined on semidi- 
rect products. We present the abstract theory, but of course historically this grew 
out of the examples, especially the heavy top and compressible flow. 

Generalities on Semidirect Products. We begin by recalling some definitions 
and properties of semidirect products. Let V be a vector space and assume that 
the Lie group G acts on the left by linear maps on V (and hence G also acts on on 
the left on its dual space V*). As sets, the semidirect product S = G©V is the 
Cartesian product S = G x V whose group multiplication is given by 

{gi,v 1 )(g 2 ,v 2 ) = {gi92,vi +giv 2 ), (2.1) 

where the action of g G G on v G V is denoted simply as gv. The identity element 
is (e, 0) where e is the identity in G. We record for convenience the inverse of an 
element: 

M-^ig-^-g^v). (2.2) 

The Lie algebra of S is the semidirect product Lie algebra, s = 0©^, whose 
bracket has the expression 

[(Zi,vi), (6,«2)] = - &t>i) , (2.3) 

where we denote the induced action of on V by concatenation, as in £ii>2. 

Below we will need the formulae for the adjoint and the coadjoint actions for 
semidirect products. We denote these and other actions by simple concatenation; 
so they are expressed as (see, e.g., Marsden, Ratiu and Weinstein [1984a,b]) 

(g,v)(S,u) = (gS,gu-(gOv), (2.4) 

and 

{g, v)(p, a) = (gp + p* v (ga),ga), (2.5) 

where (g, v) G S = G x V, (C, it) G 5 = x V, (n,a) G s* = 0* x V*, g£ = Ad g £, 
gfi = Ad*_i/i, ga denotes the induced left action of g on a (the left action of G on 
V induces a left action of G on V* — the inverse of the transpose of the action on 
V), p v '■ — ► V is the linear map given by p v (£) = £v, and p* : V* —>■ q* is its dual. 

Important Notation. For a G V*, we shall write, for notational convenience, 

p* v a = v o a G 0* , 

which is a bilinear operation in v and a. Using this notation, the above formula for 
the coadjoint action reads 

(g, v)(p, a) = (gp + vo (ga),ga). 
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We shall also denote actions of groups and Lie algebras by simple concatenation. 
For example, the g-action on g* and V* , which is defined as minus the dual map 
of the g-action on g and V respectively, is denoted by £/z and £a for £ £ g, \i £ g*, 
and a £ V * . 

Using this concatenation notation for Lie algebra actions provides the following 
alternative expression of the definition of v o a £ g*: For all v £ V, a £ V* and 
77 £ g, we define 

(77a, v) = — (v o a , 77) . 



Left Versus Right. When working with various models of continuum mechanics 
and plasmas it is convenient to work with right representations of G on the vector 
space V (as in, for example, Holm, Marsden and Ratiu [1986]). We shall denote 
the semidirect product by the same symbol S = G (§)V , the action of G on V being 
denoted by vg. The formulae change under these conventions as follows. Group 
multiplication (the analog of (|2.1|)) is given by 



(51 > «l) (52,^2) = (5132, v 2 + vig 2 ), (2.6) 
and the Lie algebra bracket on s = g(|)U (the analog of (|2.3| )) has the expression 

[(£i,«i),(6,i>2)] = ([&,&], ui6-«2£i), (2.7) 

where we denote the induced action of g on V by concatenation, as in v±^2- The 
adjoint and coadjoint actions have the formulae (analogs of (2.4) and (2.5)) 

(g,v)(£,u) = (g^iu + vOg- 1 ), (2.8) 



{g,v){jjt,a) = (jgn+(vg l )o{ag ^),ag x ), (2.9) 

where, as usual, g£ = Ad 9 £, g/s = Ad*_i//, ag denotes the inverse of the dual 
isomorphism defined by g £ G (so that g 1— > ag is a 7%M action). Note that the 
adjoint and coadjoint actions are left actions. In this case, the g-actions on g* and 
V* are defined as before to be minus the dual map given by the g-actions on g and 
V and are denoted by (because it is a left action) and a£ (because it is a right 
action) respectively. 



Lie-Poisson Brackets and Hamiltonian Vector Fields. For a left represen- 
tation of G on V the ± Lie-Poisson bracket of two functions /, k : 5* — > R is given 
by 



±(7X 



5/ 5A; 



' 5a 5fi 5a 



(2.10) 



where 5f /5/jl £ g, and 5f /5a £ V are the functional derivatives of /. The Hamilto- 
nian vector field of h : s* — * R has the expression 



v / \ /" 1* 5h \ 

X h {ji, a) = t I ad tffc/ ^ ^-j^ 00 '-^ ) 



(2.11) 
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Thus, Hamilton's equations on the dual of a semidirect product are given by 



Sh 

A = =F ad* Sh/Stl n± — oo, (2.12) 
Sh 

a = ± — a, (2.13) 

Of! 

where overdot denotes time derivative. For right representations of G on V the 
above formulae change to: 



{f,k}±(fi,a) = ±(fi 



5J_ 6k 

(fyt' 5jJL 



. 6k 6f 6f 6k . , ( | ^ . 

^ ' 5a Sfi 5a 5n 



X h (p, a) = T ^ad; A/i/1 // + ^ o a, a ^ , (2.15) 

Sh 

d = fa-. (2.17) 

Symplectic Actions by Semidirect Products. To avoid a proliferation of 
signs, in this section we consider all semidirect products to come from a left rep- 
resentation. Of course if the representation is from the right, there are similar 
formulae. 

We consider a symplectic action of 5 on a symplectic manifold P and assume 
that this action has an equivariant momentum map : P — > s*. Since V is a 
(normal) subgroup of S, it also acts on P and has a momentum map Jy : P — > V* 
given by 



where iy : V — > 5 is the inclusion v i— ► (0, u) and iy : s* — > V* is its dual. We think 
of this merely as saying that Jy is the second component of Jg. 

We can regard G as a subgroup of 5 by g *— > (g, 0). Thus, G also has a momentum 
map that is the first component of J5 but this will play a secondary role in what 
follows. On the other hand, equivariance of J5 under G implies the following relation 
for Jy: 

J v (gz) = gj v (z) (2.18) 

where we denote the appropriate action of g E G on an element by concatenation, 
as before. To prove ( |2.18j ), one uses the fact that for the coadjoint action of S on 
5* the second component is just the dual of the given action of G on V. 
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The Classical Semidirect Product Reduction Theorem. In a number of 
interesting applications such as compressible fluids, the heavy top, MHD, etc., one 
has two symmetry groups that do not commute and thus the commuting reduction 
by stages theorem of Marsden and Weinstein [1974] does not apply In this more 
general situation, it matters in what order one performs the reduction, which occurs, 
in particular for semidirect products. The main result covering the case of semidirect 
products has a complicated history with important early contributions by many 
authors, as we have listed in the introduction. The final version of the theorem as 
we shall use it, is due to Marsden, Ratiu and Weinstein [1984a,b]. 

The semidirect product reduction theorem states, roughly speaking, that for the 
semidirect product S = G (S) V where G is a group acting on a vector space V and S 
is the semidirect product, one can first reduce T*S by V and then by G and thereby 
obtain the same result as when reducing by S. As above, we let s = g (§) V denote 
the Lie algebra of S. The precise statement is as follows. 

Theorem 2.1 (Semidirect Product Reduction Theorem.) Let S = G(S)V, 
choose a = (//, a) £ Q* x V* , and reduce T*S by the action of S at a giving the coad- 
joint orbit O a through a £ 5*. There is a symplectic diffeomorphism between O a 
and the reduced space obtained by reducing T*G by the subgroup G a (the isotropy of 
G for its action on V* at the point a £ V* ) at the point fi\g a where g a is the Lie 
algebra of G a - 

Reduction by Stages. This result is a special case of a theorem on reduction 
by stages for semidirect products acting on a symplectic manifold (see Marsden, 
Misiolek, Perlmutter and Ratiu [1997] for this and more general results and see 
Leonard and Marsden [1997] for an application to underwater vehicle dynamics). 

As above, consider a symplectic action of S on a symplectic manifold P and 
assume that this action has an equivariant momentum map Js : P — > s*. As we 
have explained, the momentum map for the action of V is the map Jy : P — > V* 
given by Jy = iy o Js 

We carry out the reduction of P by S at a regular value a = (/x, a) of the 
momentum map Js for S in two stages using the following procedure. First, reduce 
P by V at the value a (assume it to be a regular value) to get the reduced space 
Pa = Jy 1 ( a )/^- Second, form the group G a consisting of elements of G that leave 
the point a fixed using the action of G on V*. One shows (and this step is not 
trivial) that the group G a acts on P a and has an induced equivariant momentum 
map J a ■ P a — ^ 0a > where Q a is the Lie algebra of G a , so one can reduce P a at the 
point n a := fi\g a to get the reduced space (P )/x„ = 3 a 1 {^a) / {G a )^ a - 

Theorem 2.2 (Reduction by Stages for Semidirect Products.) The reduced 
space (P a )iJ,a * s symplectically diffeomorphic to the reduced space P a obtained by re- 
ducing P by S at the point a = (/i, a). 

Combined with the cotangent bundle reduction theorem (see Abraham and Mars- 
den [1978] and Marsden [1992] for an exposition and references), the semidirect 
product reduction theorem is a useful tool. For example, using these tools, one 
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sees readily that the generic coadjoint orbits for the Euclidean group are cotangent 
bundles of spheres with the associated coadjoint orbit symplectic structure given by 
the canonical structure plus a magnetic term. 



Semidirect Product Reduction of Dynamics. There is a technique for re- 
ducing dynamics that is associated with the geometry of the semidirect product 
reduction theorem. One proceeds as follows: 

• We start with a Hamiltonian H ao on T*G that depends parametrically on a 
variable oo G V*. 

• The Hamiltonian, regarded as a map H : T*G x V* — > 1R is assumed to be 
invariant on T*G under the action of G on T*G x V*. 

• One shows that this condition is equivalent to the invariance of the function 
H defined on T*S = T*G x V x V* extended to be constant in the variable V 
under the action of the semidirect product. 

• By the semidirect product reduction theorem, the dynamics of H aQ reduced 
by G ao , the isotropy group of ao, is symplectically equivalent to Lie-Poisson 
dynamics on s* = q* x V*. 



• This Lie-Poisson dynamics is given by the equations ( 2.12| ) and ( |2.13| ) for the 
function h(/i,a) = H(a g , g~ x a) where \i = g~ x a g . 

3 Lagrangian Semidirect Product Theory 

Despite all the activity in the Hamiltonian theory of semidirect products, little 
attention has been paid to the corresponding Lagrangian side. Now that Lagrangian 
reduction is maturing (see Marsden and Scheurle [1993a,b]), it is appropriate to 
consider the corresponding Lagrangian question. We shall formulate four versions, 
depending on the nature of the actions and invariance properties of the Lagrangian. 
(Two of them are relegated to the appendix.) 

It should be noted that none of the theorems below require that the Lagrangian 
be nondegenerate. The subsequent theory is entirely based on variational principles 
with symmetry and is not dependent on any previous Hamiltonian formulation. We 
shall, however, show that this purely Lagrangian formulation is equivalent to the 
Hamiltonian formulation on duals of semidirect products, provided an appropriately 
defined Legendre transformation happens to be a diffeomorphism. 

The theorems that follow are modelled after the reduction theorem for the basic 
Euler-Poincare equations given earlier. However, as we shall explain, they are not 
literally special cases of it. To distinguish the two types of results, we shall use 
phrases like basic Euler-Poincare equations for the equations ( |1.57 ) and simply the 



Euler-Poincare equations or the Euler-Poincare equations with advection or the 
Euler-Poincare equations with advected parameters, for the equations that follow. 

The main difference between the left (right) invariant Lagrangians considered in 
the theorem above and the ones we shall work with below is that L and I depend 
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in addition on another parameter a £ V*, where V is a representation space for 
the Lie group G and L has an invariance property relative to both arguments. As 
we shall see below, the resulting Euler Poincare equations are not the Euler- 
Poincare equations for the semidirect product Lie algebra g (|) V* or on g (s) V, for 
that matter. 

Upcoming Examples. As we shall see in the examples, the parameter a £ V* 
acquires dynamical meaning under Lagrangian reduction. For the heavy top, the 
parameter is the unit vector in the direction of gravity, which becomes a dynami- 
cal variable in the body representation. For compressible fluids, the parameter is 
the density of the fluid in the reference configuration, which becomes a dynamical 
variable (satisfying the continuity equation) in the spatial representation. 

Left Representation and Left Invariant Lagrangian. We begin with the 
following ingredients: 

• There is a left representation of Lie group G on the vector space V and G acts 
in the natural way on the left on TG x V*: h(v g , a) = (hv g , ha). 

• Assume that the function L : TG x V* — > R is left G-invariant. 

• In particular, if ao £ V*, define the Lagrangian L ao : TG — > R by L ao (v g ) = 
L(v g ,ao). Then L ao is left invariant under the lift to TG of the left action of 
G aQ on G, where G ao is the isotropy group of ao . 

• Left G-invariance of L permits us to define I : x V* — > R by 

Kg^Vg^^ao) = L(v g ,a ). 

Conversely, this relation defines for any I : q x V* — > R a left G-invariant 
function L : TG xP-tl. 

• For a curve g(t) G G, let 

and define the curve a(t) as the unique solution of the following linear differ- 
ential equation with time dependent coefficients 

a(t) = -£(t)a(t), 

with initial condition a(0) = ao- The solution can be written as a(t) = 
g(t)~ l a Q . 

Theorem 3.1 With the preceding notation, the following are equivalent: 
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i With do held fixed, Hamilton's variational principle 

S I* L ao (g(t),g(t))dt = (3.1) 

holds, for variations 5g(t) of g(t) vanishing at the endpoints. 

ii g(t) satisfies the Euler-Lagrange equations for L ao on G. 

iii The constrained variational principle^ 

5 ( 2 l{i{t),a{t))dt = (3.2) 

Jtl 

holds on g x V* , using variations of € and a of the form 

£f = J7+ [£,"?], 5a = -rja, (3.3) 
where rj(t) £ q vanishes at the endpoints. 
iv The Euler Poincare equation^hold on q x V* 

d 51 51 51 

T - = ad? — + — oa. (3.4) 
dt d£ ? 5£ 5a y 1 



Proof. The equivalence of i and ii holds for any configuration manifold and so, 
in particular, it holds in this case. 

Next we show the equivalence of iii and iv. Indeed, using the definitions, inte- 
grating by parts, and taking into account that r\{t\) = ^(£2) = 0> we compute the 
variation of the integral to be 

f t2 I d fSl\ „ 51 51 

5 As with the basic Euler-Poincare equations, this is not strictly a variational principle in the 
same sense as the standard Hamilton's principle. It is more of a Lagrange d'Alembert principle, 
because we impose the stated constraints on the variations allowed. 

6 Note that these equations are not the basic Euler-Poincare equations because we are not 
regarding jj x V* as a Lie algebra. Rather these equations are thought of as a generalization of the 
classical Euler-Poisson equations for a heavy top, written in body angular velocity variables, as we 
shall see in the examples. Some authors may prefer the term Euler-Poisson-Poincare equations for 
these equations. 
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and so the result follows. 

Finally we show that i and iii are equivalent. First note that the G-invariance 
of L : TG x V* — ► M and the definition of a(t) = g(t)~ 1 ao imply that the integrands 



in (3A) and ( |3.2j ) are equal. However, all variations 5 git) £ TG of g(t) with fixed 
endpoints induce and are induced by variations 5£(t) G q of of the form 5£ = 
V + [£> ^1 with r/(i) S g vanishing at the endpoints; the relation between 5g(t) and 
r)(t) is given by rj(t) = git)~ 1 5git). This is the content of the following lemma 
proved in Bloch et al. [1996]. [] 

Lemma 3.2 Let g : U Cl 2 — > G 6e a smooth map and denote its partial derivatives 
by 

£(t,e) =TL g(ti£) -x(dg(t,e)/dt) 

and 

r,(t,e) =TL g(tj£) - 1 (dg(t,£)/d£). 

Then 

(3-5) 

Conversely, if U is simply connected and : U — > g are smooth functions sat- 
isfying (jS.dj) then there exists a smooth function g : U —> G such that e) = 
TL g{t:£) -i(dg(t,e)/dt) andr](t,e) = TL g[t ^ e) -i (dg(t, e)/de). 

Thus, if i holds, we define n(t) = g(t)~ 1 5g(t) for a variation <5g(i) with fixed 
endpoints. Then if we let <5£ = o(i) -1 <?(^)> we have by the above proposition <5£ = f]+ 
[i,rj\. In addition, the variation of a{t) = g(t) clq is <5a(t) = —rj(t)a(t). Conversely, 
if 6£ = f] + [£, ?y] with 77(f) vanishing at the endpoints, we define 6g(t) = g(t)r/(t) 
and the above proposition guarantees then that this Sg(t) is the general variation of 
g{t) vanishing at the endpoints. From Sa(t) = —n(t)a(t) it follows that the variation 
of g(t)a(t) = ao vanishes, which is consistent with the dependence of L ao only on 

g(t),g(t). ■ 

Cautionary Remarks. Let us explicitly show that these Euler-F 'oincare equa- 



tions (3.4) are not the Euler-P oincare equations for the semidirect product Lie al- 



gebra Q®V*. Indeed, by (|1.57| ) the basic Euler-Poincare equations 



d M 1* &l //- \ ^t^* 



dt6(Z,a) 6(£,a 

for I : q (D V* — > R become 

d 0/ „ 51 51 d 51 51 

= ad f 1 on, = — £ — , 

dt o£ ? 5£ 5a dt 5a 5a 

which is a different system from that given by the Euler-Poincare equation ( |3.4D 
and a = — £a, even though the first equations of both systems are identical. 

7 This lemma is simple for matrix groups, as in Marsden and Ratiu [1994] , but it is less elementary 
for general Lie groups. 
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The Legendre Transformation. As we explained earlier, one normally thinks of 
passing from Euler-Poincare equations on a Lie algebra q to Lie-Poisson equations 
on the dual q* by means of the Legendre transformation. In our case, we start 
with a Lagrangian on q x V* and perform a partial Legendre transformation in the 
variable £ only, by writing 



51 u \ 
M = T7, h(fi,a) 



Since 



— - £ + 

6[A 



5£ 



(3.6) 



/a 



and 5h/5a = —51 /5a, we see that ( p^|) and d(i) = - 
the minus Lie-Poisson bracket (that is, the sign + in 



-£(t)a(t) imply ( ^Tl| ) for 
pTTID ). If this Legendre 



transformation is invertible, then we can also pass from the the minus Lie-Poisson 



equations (2.11) to the Euler-Poincare equations (3.4) together with the equations 
a(t) = -£(t)a(t). 



Right Representation and Right Invariant Lagrangian. There are four ver- 
sions of the preceding theorem, the given left-left version, a left-right, a right-left 
and a right-right version. For us, the most important ones are the left-left and the 
right-right versions. We state the remaining two in the appendix. 
Here we make the following assumptions: 

• There is a right representation of Lie group G on the vector space V and G 
acts in the natural way on the right on TG x V*: (v g ,a)h = (v g h,ah). 

• Assume that the function L : TG x V* — > M is right G-invariant. 

• In particular, if ao 6 V*, define the Lagrangian L ao : TG — ► R by L aQ (v g ) = 
L(v g , ao)- Then L ao is right invariant under the lift to TG of the right action 
of G ao on G, where G ao is the isotropy group of ao- 

• Right G-invariance of L permits us to define I : q x V* — > R by 

l(y g g~ l ,aog~ l ) = L(v g ,a ). 

Conversely, this relation defines for any I : g x V* — > R a right G-invariant 
function L : TG xV'^i 

• For a curve g(t) € G, let £(t) := 5(t)g((t) _1 and define the curve a(t) as 
the unique solution of the linear differential equation with time dependent 
coefficients a(t) = —a(t)£(t) with initial condition a(0) = ao- The solution can 
be written as a(t) = aogft)^ 1 . 



Theorem 3.3 The following are equivalent: 
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i Hamilton's variational principle 



t 3 

L ao (g(t),g{t))dt = (3.7) 



holds, for variations 5g(t) of g(t) vanishing at the endpoints. 

ii g(t) satisfies the Euler-Lagrange equations for L ao on G. 

iii The constrained variational principle 

ft 2 



l(C(t),a(t))dt = (3.8) 

holds on q x V* , using variations of the form 

^ = V-[£,V], 8a = -art, ( 3 -9) 
where rj(t) £ Q vanishes at the endpoints. 
iv The Euler-Poincare equations hold on q x V* 

d 81 ,„ 51 81 . 

sr-^r*"' (3 - 10) 



The same partial Legendre transformation ( |3.6| ) as before maps the Euler- 
Poincare equations ( 3.10| ), together with the equations a = — a£ for a to the plus 
Lie-Poisson equations (2.16) and ( |2.17| ) (that is, one chooses the overall minus sign 
in these equations). 



Generalizations. The Euler-Poincare equations are a special case of the reduced 
Euler-Lagrange equations (see Marsden and Scheurle [1993b] and Cendra, Marsden 
and Ratiu [1997]). This is shown explicitly in Cendra, Holm, Marsden and Ratiu 
[1997]. There is, however, an easy generalization that is needed in some of the 
examples we will consider. Namely, if L : TG x V* x TQ and if the group G acts 
in a trivial way on TQ, then one can carry out the reduction in the same way as 
above, carrying along the Euler-Lagrange equations for the factor Q at each step. 
The resulting reduced equations then are the Euler-Poincare equations above for the 
q factor, together the Euler-Lagrange equations for the q £ Q factor. The system 
is coupled through the dependence of L on all variables. (For a full statement, 
see Cendra, Holm, Hoyle and Marsden [1997], who use this extension to treat the 
Euler-Poincare formulation of the Maxwell- Vlasov equations for plasma physics.) 



4 The Kelvin-Noether Theorem 

In this section, we explain a version of the Noether theorem that holds for solutions 
of the Euler-Poincare equations. Our formulation is motivated and designed for 
ideal continuum theories (and hence the name Kelvin-Noether), but it may also of 
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interest for finite dimensional mechanical systems. Of course it is well known (going 
back at least to the pioneering work of Arnold [1966a]) that the Kelvin circulation 
theorem for ideal flow is closely related to the Noether theorem applied to continua 
using the particle relabelling symmetry group. 

There is a version of the theorem that holds for each of the choices of conventions, 
but we shall pick the left-left conventions to illustrate the result. 

The Kelvin-Noether Quantity. We start with a Lagrangian L ao depending on a 
parameter ao € V* as above. We introduce a manifold C on which G acts (we assume 
this is also a left action) and suppose we have an equivariant map K, : C x V* — > g** . 

As we shall see, in the case of continuum theories, the space C will be a loop 
space and ()C(c,a),u) for c G C and /i £ g* will be a circulation. This class of 
examples also shows why we do not want to identify the double dual g** with g. 

Define the Kelvin-Noether quantity I : C x g x V* — > IR by 



We are now ready to state the main theorem of this section. 

Theorem 4.1 (Kelvin-Noether.) Fixing Co 6 C, let £(t),a(t) satisfy the Euler- 
Poincare equations and define g(t) to be the solution of g(t) = g(t)£(t) and, say, 
g(0) = e. Let c(t) = g{t)~ l c a and I(t) = I(c(t),£(t),a(t)). Then 



Proof. First of all, write a(t) = g(t) 1 ao as we did previously and use equivariance 
to write I(t) as follows: 



The g^ 1 pulls over to the right side as g (and not g~ l ) because of our conventions 
of always using left representations. We now differentiate the right hand side of this 
equation. To do so, we use the following well known formula for differentiating the 
coadjoint action (see Marsden and Ratiu [1994], page 276): 




(4.1) 




(4.2) 




-[g(t)u(t)} = g(t) - ad* (t) u(t) + -u(t) 



where, as usual, 
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Using this coadjoint action formula and the Euler-Poincare equations, we obtain 



dt 



d_ 

di 



K(cQ,a ),g(t) 
d 



5l_ 

5£ 



(£(t),a(t)) 



K(co,ao),g(t) 
K(co,ao),g(t) 
5(t) _1 /C(c ,a 
K(c(t),a(t)), 



(t) 



m,a(t)) 
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51 51 



51 

— o a 
5a 



51 

— o a 
5a 



51 

— o a 
5a 



where, in the last steps, we used the definitions of the coadjoint action, as well as 
the Euler-Poincare equation (pTj) and the equivariance of the map /C. ■ 



Corollary 4.2 For the basic Euler-Poincare equations, the Kelvin quantity I(t), 
defined the same way as above but with I : C x q — > R, is conserved. 

For a review of the standard Noether theorem results for energy and momentum 
conservation in the context of the general theory, see, e.g., Marsden and Ratiu [1994]. 



5 The Heavy Top 

In this section we shall use Theorem |3.1| to derive the classical Euler-Poisson equa- 
tions for the heavy top. Our purpose is merely to illustrate the theorem with a 
concrete example. 



The Heavy Top Lagrangian. The heavy top kinetic energy is given by the 
left invariant metric on 50(3) whose value at the identity is (fti,^) = I^i • 
where fii,02 £ H^ 3 are thought of as elements of so(3), the Lie algebra of SO(3), 
via the isomorphism ft £ R 3 i— > ft £ so (3), fiv := fl x v, and where I is the 
(time independent) moment of inertia tensor in body coordinates, usually taken as 
a diagonal matrix by choosing the body coordinate system to be a principal axes 
body frame. This kinetic energy is thus left invariant under the full group £0(3). 
The potential energy is given by the work done in lifting the weight of the body to the 
height of its center of mass, with the direction of gravity pointing downwards. If M 
denotes the total mass of the top, g the magnitude of the gravitational acceleration, 
X the unit vector of the oriented line segment pointing from the fixed point about 
which the top rotates (the origin of a spatial coordinate system) to the center of mass 
of the body, and t its length, then the potential energy is given by — MgtTtT^ez ■ x, 
where e3 is the axis of the spatial coordinate system parallel to the direction of 
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gravity but pointing upwards. This potential energy breaks the full SO (3) symmetry 
and is invariant only under the rotations S 1 about the e3~axis. 

However, for the application of Theorem |3.1| we are supposed to think of the 
Lagrangian of the heavy top as a function on TSO(3) x K 3 -> 1. That is, we need 
to think of the potential energy as a function of (%, v) 6 TSO(3) x R 3 . This means 
that we need to replace the vector giving the direction of gravity e3 by an arbitrary 
vector v £ R 3 , so that the potential equals 

i7(ur,v) = M 9 m-\- X - 

Thought of this way, the potential is 50(3)-invariant. Indeed, if R' G SO{3) is 
arbitrary, then 

[/(R'u R ,R'v) = Mgt (R'R) x R'v x 
= MglTT l v-x 
= ?7(ur,v) 



and the hypotheses of Theorem 3.1 are satisfied. Thus, the heavy top equations of 
motion in the body representation are given by the Euler-Poincare equations ( |3.4[ ) 
for the Lagrangian I : so (3) x R 3 -> K. 

The Reduced Lagrangian. To compute the explicit expression of I, denote by 
fi the angular velocity and by II = Ifi the angular momentum in the body rep- 
resentation. Let r = R x v; if v = e3, the unit vector pointing upwards on the 
vertical spatial axis, then T is this unit vector viewed by an observer moving with 
the body. The Lagrangian I : so (3) x R 3 — > R is thus given by 

l(fl,T) = L(R" 1 n R ,R" 1 v) 

= in-n + [/(R- 1 u R ,R- 1 v) 
= hi-n + Mgiv-x- 

The Euler Poincare Equations. It is now straightforward to compute the 
Euler-Poincare equations. First note that 

Since 

ad£j n = nxfi, vor = -Ixv, 

and 

nr = -r x n, 

the Euler-Poincare equations are 

ri = n x n + MglT x x, 
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which are coupled to the T evolution 

r=rxfi. 

This system of two vector equations comprises the classical Euler-Poisson equations, 
which describe the motion of the heavy top in the body representation. 



The Kelvin-Noether theorem Let C = Q and let K : C x V* -> g** = q be the 

map (W, r) i — > W. Then the Kelvin-Noether theorem gives the statement 

j t (W,U) = Mg£ (W,rxx) 

where W(i) = R(i) _1 w; in other words, W(i) is the body representation of a 
space fixed vector. This statement is easily verified directly. Also, note that 
(W, II) = (w, 7r), with 7r = R(t)II, so the Kelvin-Noether theorem may be viewed 
as a statement about the rate of change of the momentum map of the system (the 
spatial angular momentum) relative to the full group of rotations, not just those 
about the vertical axis. 



6 The Euler-Poincare Equations 
in Continuum Mechanics 

In this section we will apply the Euler-Poincare equations in the case of continuum 
mechanical systems. We let D be a bounded domain in M. n with smooth boundary 
&D (or, more generally, a smooth compact manifold with boundary and given vol- 
ume form or density). We let Diff (T>) denote the diffeomorphism group of T> of some 
given Sobolev class. If the domain T> is not compact, then various decay hypotheses 
at infinity need to be imposed. Under such conditions, Diff(2?) is a smooth infi- 
nite dimensional manifold and a topological group relative to the induced manifold 
topology. Right translation is smooth but left translation and inversion are only 
continuous. Thus, Diff (2?) is not actually a Lie group and the previous theory does 
not apply, strictly speaking. Nevertheless, if one uses right translations and right 



representations, the Euler-Poincare equations of Theorem 3J3 do make sense, as a 
simple verification shows. We shall illustrate below such computations, by verifying 
several key facts in the proof. 

Let X(T>) denote the space of vector fields on T> of some fixed differentiability 
class. Formally, this is the right Lie algebra of Diff (T>) , that is, its standard left 
Lie algebra bracket is minus the usual Lie bracket for vector fields. To distinguish 
between these brackets, we shall reserve in what follows the notation [u, v] for the 
standard Jacobi-Lie bracket of the vector fields u, v £ 3C(T>) whereas the notation 
ad u v := — [u, v] denotes the adjoint action of the left Lie algebra on itself. 

We also let X(T>)* denote the geometric dual space of X(T>), that is, X(V)* := 
il 1 (D) ®Den(P), the space of one-form densities on T>. If a (%>m 6 f2 (£>) (g>Den(2?), 
the pairing of a (g> m with u € X(T>) is given by 



(a<S>m, u) = / a ■ um (6-1) 



v 
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where a ■ u is the standard contraction of a one-form with a vector field. For 
u G X{V) and a m G X(D)*, the dual of the adjoint representation is defined by 



to 



(ad* (a <g) m), v) = — / a-[u,v] 
Jv 

and its expression is 

ad* (a ® m) = (i? u a + (div m u)a) ®m = £ u (a ® m) , (6-2) 

where div m u is the divergence of u relative to the measure to, that is, £ u m = 
(div m u)m. Hence if u = u l d/dx l , a = aidx 1 , the one-form factor in the preceding 
formula for ad* (a <g) m) has the coordinate expression 

■ dot ■ 3n? \ ■ id d*uP \ 

u3 l^ + a il^ + ( div - U H Uc* = + a j7 r-) dx l , (6.3) 



the last equality assuming that the divergence is taken relative to the standard 
measure m = d n x in W 1 . (On a Riemannian manifold the metric divergence needs 
to be used.) 

Throughout the rest of the paper we shall use the following conventions and 
terminology for the standard quantities in continuum mechanics. Elements of T> 
representing the material particles of the system are denoted by X; their coordinates 
X A ,A = 1, ...,n may thus be regarded as the particle labels. A configuration, 
which we typically denote by i], is an element of Diff (T>). A motion r/t is a path in 
Diff(P). The Lagrangian or material velocity \~(X,t) of the continuum along 
the motion rjt is defined by taking the time derivative of the motion keeping the 
particle labels (the reference particles) X fixed: 



Vt(X), 

x 



the second equality being a convenient shorthand notation of the time derivative for 
fixed X. 

Consistent with this definition of velocity, the tangent space to Diff(P) at rj G 
Diff(P) is given by 

T v Diff (V) = {V ri : V ->■ TV \ V n (X) G T v{x) V}. 

Elements of T T] Diff (I?) are usually thought of as vector fields on V covering r\. The 
tangent lift of right and left translations on TDiff(£>) by ip G Diff (V) have the 
expressions 

V v (p :=T r) R lf (V ri ) = V v o<p and ipV^ := T^L^iy^) = Tip o W, . 

During a motion r/t, the particle labeled by X describes a path in T> whose 
points x(X, t) := r/ t (X) are called the Eulerian or spatial points of this path. 
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The derivative v(x, t) of this path, keeping the Eulerian point x fixed, is called the 
Eulerian or spatial velocity of the system: 



v(x,t):=V(X,t):= - Vt (X). 



Thus the Eulerian velocity v is a time dependent vector field on T>: v< £ 
where v t (x) := v(x,i). We also have the fundamental relationship 



The representation space V* of Diff (V) in continuum mechanics is often some 
subspace of T(2?) <g> Den(P), the tensor field densities on T> and the representation 
is given by pull back. It is thus a right representation of Diff (P) on Z(T>) (g)Den(X>). 
The right action of the Lie algebra X(V) on V* is given by av := £ v a, the Lie 
derivative of the tensor field density a along the vector field v. 

The Lagrangian of a continuum mechanical system is a function L : T Diff (I?) x 
V* — > R which is right invariant relative to the tangent lift of right translation of 
Diff (I?) on itself and pull back on the tensor field densities. 

Thus, the Lagrangian L induces a function I : X(V) xP-»l given by 



where v 6 3C(T>) and a £ V* C T(2?) <8> Den(P), and where 77* a denotes the pull 
back of a by the diffeomorphism 77 and v is the Eulerian velocity. The evolution of 
a is given by the equation 



The solution of this equation, given the initial condition ao, is a(t) = (pt*ao, where 
the lower star denotes the push forward operation and ip t is the flow of v. 

Advected Eulerian quantities are defined in continuum mechanics to be those 
variables which are Lie transported by the flow of the Eulerian velocity field. Using 
this standard terminology, the above equation states that the tensor field density a 
(which may include mass density and other Eulerian quantities) is advected. 

As remarked, V* C T(£>) <8> Den(D). On a general manifold, tensors of a given 
type have natural duals. For example, symmetric covariant tensors are dual to 
symmetric contravariant tensor densities, the pairing being given by the integration 
of the natural contraction of these tensors. Likewise, fc-forms are naturally dual to 
(n — fc)-forms, the pairing being given by taking the integral of their wedge product. 

The operation o between elements of V and V* producing an element of X(V)* 
introduced in section 2 becomes 



Jv 

where v ■ £ u a denotes the contraction, as described above, of elements of V and 
elements of V*. (These operations do not depend on a Riemannian structure.) 



V t = v t o rit , 



where V t {X) :=V(X,t). 



l(v, a) = L(v o 77, rfa) 



a = — £ v a. 




(6.4) 



34 



For a path r/t S DifF(P) let v(x,t) be its Eulerian velocity and consider as in 
the hypotheses of Theorem |3.3j the curve a(t) with initial condition ao given by the 
equation 

a + £ v a = 0. (6-5) 



Let L ao (V) := L(V, ao). We can now state Theorem in this particular, but very 
useful, setting. 

Theorem 6.1 (Euler Poincare Theorem for Continua.) Consider a path r\t 
in Diff(D) with Lagrangian velocity V and Eulerian velocity v. The following are 
equivalent: 

i Hamilton's variational principle 



5 I 2 L(X,V t (X),ao(X)) dt = 



(6.6) 



holds, for variations Srjt vanishing at the endpoints. 

ii rjt satisfies the Euler-Lagrange equations for L ao on Diff(D).[] 

iii The constrained variational principle in Eulerian coordinates 

ft 2 



' / l(y,a) 
Jti 



dt = (6.7) 



holds on 3L(T>) x V* , using variations of the form 



du 

5v = — + [v,u], 5a = -£ u a, (6.8) 



where ut = 5rjt oi| t 1 vanishes at the endpoints. 
iv The Euler~P oincare equations for continua 



8 51 ,„ 51 51 n 51 51 

di^ = -^^ + Ta «* = -^ 1 - + Ta «a, (6.9) 



hold, where the o operation given by ( \6.Q ) needs to be determined on a case by 
case basis, depending on the nature of the tensor a. (Remember that 51 /5v is 
a one-form density.) 



We do not write these equations explicitly since to do so would require either a coordinatization 
of the diffeomorphism group, which is not easy to give explictly, or requires more structure, such 
as an affine connection on this group. Certainly writing the equations formally, imagining that rj 
and 77 form valid coordinates in which the Euler-Lagrange equations hold is not correct. 
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Remarks. 



1. Of course, this theorem can be proven directly by imitating the proof of The- 



orem 3.1 with appropriate modifications for right representations and right 
actions. For those used to the more concrete language of continuum mechan- 
ics as opposed to that of Lie algebras, the following string of equalities shows 
that iii is equivalent to iv: 







5 I l(v,a)dt 
hi 



t-2 



61 . 51 . . „ 
— • ov + — • da I at 
ov oa 



t-2 



51 fdu 

5^'{m- advU 



51 p ' 
oa 



u 



fi 



d_5l_ 



ad v - — h — o a 
ov 6a 



dt 
dt. 



(6.10) 



4. 



2. Similarly, one can deduce by hand the form ( p.8|) of the variations in the 
constrained variational principle ( |6.7| ) by a direct calculation. This proceeds 
as follows. One writes the basic relation between the spatial and material 
velocities, namely v(x, t) = ?)(^ 1 (x), t). One then takes the variation of 
this equation with respect to r\ and uses the definition u(x,t) = 5'i](r]^ 1 (x), t) 
together with a calculation of its time derivative. Of course, one can also do 
this calculation using the inverse map rj^~ instead of rj as the basic variable, 
see Holm, Marsden, and Ratiu [1986], Holm [1996a, b]. 

3. As we mentioned in the context of perfect fluids, the preceding sort of calcu- 
lation for oV in fluid mechanics and the interpretation of this restriction on 
the form of the variations as the so-called Lin constraints is due to Bretherton 
[1970]. 

The coordinate expressions for (51 /5a) o a required to complete the equations 
of motion are given in the next section for several choices of ao(X) in three 
dimensions. Namely, we shall discuss the choices corresponding to scalar s, 
one-forms, two-forms, densities in three dimensions, and symmetric tensors. 
In the equations of motion, all of these quantities will be advected. 

As with the general theory, variations of the action in the advected tensor 
quantities contribute to the equations of motion which follow from Hamilton's 
principle. At the level of the action I for the Euler-Poincare equations, the 
Legendre transform in the variable v alone is often nonsingular, and when it 
is, it produces the Hamiltonian formulation of Eulerian fluid motions with a 
Lie-Poisson bracket defined on the dual of the semidirect product algebra of 
vector fields acting amongst themselves by Lie bracket and on tensor fields 
and differential forms by the Lie derivative. This is a special instance of the 
more general facts for Lie algebras that were discussed earlier. 

6. As mentioned earlier, in the absence of the tensor fields a and when I is the 
kinetic energy metric, the basic Euler-Poincare equations are the geodesic 



5. 
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spray equations for geodesic motion on the diffeomorphism group with respect 
to that metric. See, e.g., Arnold [1966a], Ovsienko and Khesin [1987], Zeitlin 
and Kambe [1993], Zeitlin and Pasmanter [1994], Ono [1995a, 1995b] and 
Kouranbaeva [1997] for details in particular applications of ideal continuum 
mechanics. 

Remarks on the inverse map and the tensor fields a for fluids. In the 

case of fluids in the Lagrangian picture, the flow of the fluid is a diffeomorphism 
which takes a fluid parcel along a path from its initial position A, in a "reference 
configuration" to its current position x in the "container", i.e., in the Eulerian 
domain of flow. As we have described, this forward map is denoted by r\ : X i— > x. 
The inverse map r?^ 1 : x i— > X provides the map assigning the Lagrangian labels 
to a given spatial point. Interpreted as passive scalars, these Lagrangian labels 
are simply advected with the fluid flow, X = 0. In the Lagrangian picture, a 
tensor density in the reference configuration cto(A) (satisfying do (A) = 0) consists 
of invariant tensor functions of the initial reference positions and their differentials. 
These tensor functions are parameters of the initial fluid reference configuration 
(e.g., the initial density distribution, which is an invariant n-form). 
When viewed in the Eulerian picture as 



the time invariant tensor density ao(X) in the reference configuration acquires ad- 
vective dynamics in the Eulerian picture, namely 



where £ v denotes Lie derivative with respect to the Eulerian velocity field v(x,i). 
This relation results directly from the well known Lie derivative formula for tensor 
fields. (See, for example, Abraham, Marsden and Ratiu [1988].) 

Mapping the time invariant quantity ao (a tensor density function of X) to the 
time varying quantity at (a tensor density function of x) as explained above is a 
special case of the general way we advect quantities in V* in the general theory. 
Specifically, we can view this advection of at as being the fluid analogue of the 
advection of the unit vector along the direction of gravity (a spatially fixed quantity) 
by means of the body rotation vector in the heavy top example. 

Consistent with the fact that the heavy top is a left invariant system while 
continuum theories are right invariant, the advected tensor density at is a spatial 
quantity, while the advected direction of gravity is a body quantity. If we were to 
take the inverse map rj -1 as the basic group variable, rather than the map r], then 
continuum theories would also become left invariant. 



a t (x) := {r) t *a )(x) 



i.e., 



a (X) := (v* t a t )(X) 



(Vu'aoXX), 
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The continuity equation for the mass density. We will need to impose an 
additional assumption on our continuum theory. Namely, we assume that amongst 
the tensor densities being advected, there is a special one, namely the mass density. 
This of course is a tensor density that occurs in all continuum theories. We denote 
this density by pd n x and it is advected according to the standard principles discussed 
above. Thus, p satisfies the usual continuity equation: 

d 

—p + div(pv) = 0. 

In the Lagrangian picture we have pd n x = po(X)d n X , where po(X) is the mass 
density in the reference configuration. It will also be convenient in the continuum 
examples below to define Lagrangian mass coordinates £(X) satisfying pd n x = d n l 
with £ = 0. (When using Lagrangian mass coordinates, we shall denote the density 
p as D.) We assume that p (or D) is strictly positive. 



The Kelvin-Noether Circulation Theorem. Let C be the space of continuous 
loops 7 : S 1 — ► V in V and let the group Diff(P) act on C on the left by (77,7) 6 
Diff (T>) x C 1— »• 7/7 G C, where 777 = 77 o 7. 

Next we shall define the circulation map K : C x V* — > X(V)**. Given a one 
form density a £ X* we can form a one form (no longer a density) by dividing it by 
the mass density p; we denote the result just by a/ p. We let K, then be defined by 

(£(7, a), a) = (6.11) 

The expression in this definition is called the circulation of the one-form ot/p 
around the loop 7. 

This map is equivariant in the sense that 

{K(j) o 7, 77*0), 77*0;) = (£(7, a), a) 

for any 7/ £ Diff (D). This is proved using the definitions and the change of variables 
formula. 

Given the Lagrangian I : 3C(T>) x V* — > M., the Kelvin-Noether quantity is given 
by (|4.1|) which in this case becomes 



/(7 ' V ' fl) = /J^- 



With these definitions of K, and /, the statement of Theorem |4.1| becomes the clas- 
sical Kelvin circulation theorem. 



Theorem 6.2 (Kelvin Circulation Theorem.) Assume that v(x,t) satisfies the 
Euler-Poincare equations for continua: 

d fSl\ n {5l\ 61 
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and a satisfies the advection relation 

da 
dt 



+ £ v a = 0. 



Letrjt be the flow of the Eulerian velocity field v , that is, \ t = (drjt / dt) o rj t 1 . Define 
It ■= Vt ° 7o and I(t) := I{jt,v t , a t ). Then 



dt 



I(t) 



15l_ 

p 5a 



o a . 



In this statement, we use a subscript t to emphasise that the operations are done at 
a particular t and to avoid having to write the other arguments, as in at(x) = a(x, t); 
we omit the arguments from the notation when convenient. Due to the importance 
of this theorem we shall give here a separate proof tailored for the case of continuum 
mechanical systems. 

Proof. First we change variables in the expression for I(t): 



1_ si_ 

p t 5v 



1 * 

70 P0 



Next, we use the Lie derivative formula 

-^ {r]*a t ) = rf t (J^a t + £„at\ , 

for an arbitrary one-form density at- This formula gives 

51 



5l_ 
Sv 



d T/ , 



<5v 



- i ~rf 

dt J-yo P0 1 
I LI ( v * 

/ 70 po dt y 1 

I -rf~ 

By the Euler-Poincare equations, this becomes 



5l_ 

51 



- I — 1 - / v ( ^ 
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again by the change of variables formula. 



'51 




'51 


— oa 


= <f - 


— oa 


da 


Jit Pt 


da J 



Corollary 6.3 (Kelvin-Noether form.) Since the last expression holds for every 
loop 7 t , we may write it as 







dt 



1 51 151 



— + £ v )-— = -— oa. 



p 5v p5a 



(6.12) 



This is the Kelvin-Noether form of the Euler-Poincare equations for ideal con- 
tinuum dynamics. 
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7 Applications of the Euler— Poincare Theorem to Con- 
tinua 



Variational formulae in three dimensional Euclidean coordinates. We 

compute explicit formulae for the variations 5a in the cases that the set of tensor 
fields a consists of elements with the following coordinate functions in a Euclidean 
basis on M 3 , 

a e {b, A • dx, B • dS, D d 3 x, S ab dx a <g> dx b } . (7.1) 

These are the tensor fields that typically occur in ideal continuum dynamics. Here, 
in three dimensional vector notation, we choose B = curl A and d(A ■ dx) = B • dS. 
In Euclidean coordinates on M 3 , this is d{A k dx k ) = A k jdx^ A dx k = \tij k B l dx^ A 
dx k , where is the completely antisymmetric tensor density on R 3 with ei23 = 
+1. (The two form B-dS = d(A ■ dx) is the physically interesting special case of 
B^jdx 1 Adx k , in which B^j = A^j, so that V • B = 0.) 

We have seen that invariance of the set a in the Lagrangian picture under the 
dynamics of v implies in the Eulerian picture that (Jj + =£v) a = 0, where £ v denotes 



Lie derivative with respect to the velocity vector field v. According to Theorem 6.1 
equation ( |6.9| ) , the variations of the tensor functions a at fixed x and t are also given 
by Lie derivatives, namely da = —£ u a, or 

5b = -£ u b = -u-V6, 

JA-dx = -£ u (A-dx) = -((u- V)A + A,W) -dx 

= (u x curl A — V(u • A)) • dx , 

5B-dS = -£ u (B • dS) = (curl (u x B)) • dS = d{5A • dx) , 

5D d 3 x = -£ u (Dd 3 x) = -V • (Du) d 3 x, 

5 S ab dx a ®dx b = -£ u (Sab dx a <g) dx b ) 

= -{u k S ab , k + S kb u k a + S ka u k b )dx a ®dx b . (7.2) 
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Hence, Hamilton's principle with this dependence yields 



= 8 J dt l(v;b,A,B,D,S ab ) 



dt 
dt 



jT • 5v + — 56 + — 5D + — • 5 A + — • 5B + — - 5S ab 
5v 56 5-D 5A 5B 5S ab 



51 (d 



5v \ dt 



ad v u 



(5/ (5/ 
+ — ■ (u x curl A - V(u • A)) + — • (curl(u x B)) 
8 A 5B 



81 

5S a b 



U k S a b,k + SkbU k a + Sk a U k b 



dt 



u 



d 51 , m 61 61 , _ _ 5Z 



5/ 



+ u« 



— - x curl A + A div — - + B x curl ■— 
5A 5A 5B 

~To S a b,k 4~ (To Skb),a + ("Fo ^fca) 
0O o 6 Ob a b Ob a b 



(7.3) 



u • 

dt u • 



a 5i d * 5/ 5/ 

5t 5v v 5v 5a 



a_ \ _5/_ ' 

at / 5v 5a 



where we have consistently dropped boundary terms arising from integrations by 
parts, by invoking natural boundary conditions. Thus, for the set of tensor fields a 
in equation Q7.1| ) we have the following Euclidean components of ^oa, 



(51 



5/ n / 5/ 



5b 



+ 



5D 



61 1 A A J- 61 TD 1 5/ 

— — x curl A + A div — - + B x curl — — 
5A 5A 5B 



61 



S a bh + 



51 



Skb 



+ 



81 



'Ska 



(7.4) 



5 S a b ' \6S a b J a \5S a b 
Stress tensor formulation. For example, if we assume a Lagrangian in the form 
Z(v; b, A, B, D, S ab ) = [ d 3 x £(v, Vv, 6, A, B, D, S ab ), (7.5) 



where £ is a given function, then we may use equation ( |7.4[ ) to express the Euler- 
Poincare equations for continua Q6.9| ) in this case in the momentum conservation 
form, 



— — -ad*— — oa 
dt 5v v 5v 5a 



at dx* i ' 



(7.6) 
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with momentum density components mi, i = 1,2,3 denned by 



m .^Sl L= dC__d_(dC^\ 
% 5v l dv l dx k \ dv l k J 

and stress tensor T- given by 

rpj j dC k 9C a , dC pi 9C a 9C a 

Li ~ mV dv) V >' dAj l + dB* dS jb ib dS a j m 

Here, in the calculation of T- , we have used the coordinate expression ( fO| ) for 

ad* (a (g) m). 

Kelvin-Noether form. The Euclidean components of the Euler-Poincare equa- 
tions for ideal continua may also be summarized in Kelvin-Noether form ( |6.12 ) for 



advected tensor fields a in the set ( |7.l[) . We adopt the notational convention of the 
circulation map K. in equation ( 6.11[ ) that a one form density can be made into a one 



form (no longer a density) by dividing it by the mass density D to produce, e.g., the 
one form in Euclidean components j=j-j^dx % from the one form density J^. With a 
slight abuse of notation (but in accord with the usual physics conventions) we write 
the former coordinate expression as -^J^ • dx. We also denote the Lie-derivative 
relation for the continuity equation as (d/dt + £ w )Dd 3 x = 0. Then, the Euclidean 
components of the Euler-Poincare equations for continua in ( [f.3D are expressed in 
Kelvin-Noether form ( |6.12| ) as 

d n \ { I 51 \ 1 51 , / 51 \ , 

— + £ v ) dx) H V6 • dx - V — • dx 

dt J\D5v J D5b \5D J 

1 f 51 , A 4 SI \ , 1 /„ , SI \ , 

+ — -— x curl A - A div— — • dx - — B x curl — ) ■ dx 
D \ 6 A 5 A ) D \ 5B ) 

+ ^ ( ^S ab , k - (^S kb ), a - (^-S ka ) ib ) dx k = 0, (7.9) 
L> \5Sab ob ab db ab J 

where the components of the variational derivatives of the Lagrangian I are to be 
computed according to the usual physics conventions, i.e., as components of Frechet 
derivatives as in equation ( |1.50| ). In physical applications, the advected Eulerian ten- 



sor fields a in (7T) represent the buoyancy b (or specific entropy, for the compressible 
case), magnetic vector potential A, magnetic field intensity B, mass density D, and 
Cauchy-Green strain tensor S ab , respectively. Formula (|7.9| ) is the Kelvin-Noether 
form of the equation of motion for ideal continua in Euclidean coordinates. This 
Euclidean component formula is especially convenient for direct calculations in fluid 
dynamics, to which we turn our attention next. 
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Eulerian motion equation for an ideal incompressible fluid. In the Eulerian 
velocity representation we consider fluid motion in an n-dimensional domain and 
define the reduced action © re d and reduced Lagrangian l(v,D) by 

S red = J dtl = J dt J d n x \d\v\ 2 -p(D- 1) . (7.10) 

This action produces the following variations at fixed x and t 

1 SI 51 1 , ,9 SI , „ ,„ _ H . 



Hence, from equation ( |7.9[ ) for Hamilton principles of this type we find the Eulerian 
motion equation, 

l + ^)(^ dX )" V (^)" iX = ' 01 ^ + (vV)v + V P = 0, 

(7.12) 

for "natural" boundary conditions, n • v = on the boundary, where n is the 
boundary's outward unit normal vector. This is the Eulerian motion equation for 
an incompressible fluid in n dimensions. The constraint D = 1 (volume or mass 
preservation) is imposed by varying the Lagrange multiplier p, the presure. Incom- 
pressibility then follows from substituting D = 1 into the Lie-derivative relation for 
D, which closes the ideal incompressible fluid system, 

(9 \ 3D 

- + £ v \Dd 3 x = 0, i.e., _ = -V-(£>v). (7.13) 

This relation, together with the constraint D = 1 gives incompressibility of the flow, 
V-v = 0. 

Remark on Lagrangian mass coordinates. An alternative way to treat Hamil- 
ton's principle for the action ( 7.10| ) is to perform variations at fixed x and t of the 



inverse maps x ^ £, described by the Lagrangian mass coordinate functions ^(x, t), 
A = 1,2,... ,n, which determine v and D by the formulae (in which one sums on 
repeated indices) 

d J^ = -jDt, D? = ^-, D = det(Df). (7.14) 

As discussed above, the relation of mass coordinates t to the usual Lagrangian 
coordinates X is given by a change of variables in the fluid reference configuration 
to make po(X)d n X = d n l. Variation of an action of the form © re d(v, D) with 
respect to i A with p imposing volume preservation then yields (Holm, Marsden, and 
Ratiu [1986], Holm [1996a]), 

-dp(D-l)], (7.15) 
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where d/dt = d/dt + (v • V) is the material derivative of Eulerian quantities and we 
again invoke natural boundary conditions when integrating by parts. 

Hence, the vanishing of the coefficient of 5£ A in the variational formula (7.15) 



recovers the Euler-Poincare equation ( 7.12| ) for the Eulerian fluid velocity, v, by 



stationarity of the action (7.10) with respect to variations of the Lagrangian mass 
coordinates £ A (x, t). Similar arguments based on stationary variations of the action 
with respect to the Lagrangian mass coordinates £ A at fixed x, t will also recover 



the more general Euler-Poincare equations (|7.9|) from actions which depend on the 



velocity v and the advected quantities in equation (|7.1[) through their dependence 
on the £ A (x,t). 



Adiabatic compressible MHD. In the case of adiabatic compressible magne- 
tohydrodynamics (MHD), the action in Hamilton's principle is given by 

S red = j dtl = j dtd 3 x |v| 2 -De{D,b) - ^|B| 2 ^ , (7.16) 

where e(D, b) is the fluid's specific internal energy, whose dependence on the density 
D and specific entropy b is given as the "equation of state" and which for an isotropic 
medium satisfies the thermodynamic first law in the form de = —pd(l/D) + Tdb 
with pressure p(D, b) and temperature T(D, b). The variation of I in ( |7.16j ) is 

5& rcd = J dtd 3 x Dv-S-v- DT5b + Q|v| 2 - h\ 8D -B • 5B. (7.17) 

The quantity h = e + p/ D denotes the specific enthalpy, which thus satisfies dh = 
(1/D)dp+Tdb. The Euler-Poincare formula in the Kelvin-Noether form ( |7.9[) yields 
the MHD motion equation as 

^ + (v • dx) - Tdb + -^B x curl B-dx-d Q|v| 2 - h \ =0, (7.18) 

or, in three dimensional vector form, 

dv 11 

— + (v • V)v + — Vp + — B x curl B = 0. (7.19) 

By definition, the advected variables {b, B, D} satisfy the the following Lie-derivative 
relations which close the ideal MHD system, 

d \ db 
- + £ v j6 = 0, or ^ = -v-V6, 

d \ 9B 

— + £ v B • dS = 0, or — = curl (v x B), 
dt J dt y ' 

d . „ \ „ ;3 „ „ dD 



— + £ v j Dd d x = 0, or _ = -V-(r>v), (7.20) 

and the function p(D, b) = D 2 de/dD is specified by giving the equation of state of 
the fluid, e = e(D, b). If the condition V • B = holds initially, then it holds for all 
time; since this constraint is preserved by the dynamical equation for B. 
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Adiabatic magneto-elastodynamics. When nonlinear elasticity is also a factor 
in the MHD evolution, there is an additional Lie-derivative relation, 



^ + £v) (S ab dx a dx b ) = , (7.21) 

leading to the dynamical equation for the advected Cauchy-Green strain tensor S a b 
(which measures nonlinear strain in spatial coordinates), 

§- t S ab = ~ (v k S ab , k + S kb v k a + S ka v k h ) . (7.22) 

In this case, additional stress terms appear in the motion equation for v that arise 
from the dependence of the specific internal energy e(D,b, S a b) on the Cauchy- 
Green strain tensor S a b in the MHD action (|7.16| ) when the elasticity of the medium 
is involved. The stress tensor per unit mass a ab is determined from the equation 
of state of such an magneto-elastic medium by the Doyle-Erickson formula a ab = 
de/dS a b. The Euler-Poincare equation ( [7.9| ) for ideal magneto-elasticity is then to 
be 



r\ 11 

-jfi + + pP,i + D&iBit ~ Bij) ~ (° ab Sib),a - (v ab S ia ), b = , (7.23) 



where we have used the specific enthalpy relation for an elastic medium, dh — Tdb = 
D~ l dp + (j ab dS a b- Thus, adiabatic magneto-elastodynamics summons all of the 
advected quantities in equation ( |7.l[) and makes use of the entire Euler-Poincare 
equation ( |7.9| ). 



Adiabatic compressible Maxwell fluid dynamics via the Kaluza-Klein con- 
struction. An adiabatic Maxwell fluid (MF) with (nonrelativistic) Eulerian fluid 
velocity v, density D, specific entropy b and pressure p(D, b) satisfies the following 
system of equations, 

— + (v • V)v + -Vp(D, b) = ± (E + v x B), 
at Dm 

g-V.tf*), | = -(v.V) t , 

VE = — D, ? = --Ov|VxB, 

m at m 

V • B = 0, — = - V x E. (7.24) 
dt y ' 

This system consists of the motion equation for a charged fluid moving under the 
combined effects of pressure gradient and Lorentz forces; the continuity equation for 
the mass density D; advection of the specific entropy, b; and Maxwell's equations for 
the electromagnetic fields E and B in the moving fluid medium, whose polar izability 
and magnetization are neglected for simplicity. (For the physically more realistic 
treatment of moving media with electromagnetic induction in a similar framework, 
including relativistic effects, see Holm [1987].) The equations for D and b are the 
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familiar advection laws. The coupling constant q/m is the charge-to-mass ratio of 
the fluid particles, and the electric and magnetic fields E and B are defined in terms 
of the scalar and vector potentials $ and A by 



B=-£-V., 



B = V x A. 



(7.25) 



In the MF equations ( 7.24| ), charged fluid motion is the source for the electromag- 
netic fields which act self-consistently upon the fluid through the Lorentz force. We 
shall show that equations ( [7.24 ) are Euler-Poincare equations for the gauge invariant 
action of "Kaluza-Klein" form given by 



6 



red 



dtl 



dtd A x 



-D\w\ 2 + -D(A ■ v - $ + Co) 2 



De(D,b) 



1 

+ 2 



# A „~ 2 



-IV x A| 
2 1 1 



(7.26) 



where u = d6/dt + v • V# for a gauge field 9 and e(D,b) is the fluid's specific 
internal energy, which satisfies the first law of thermodynamics in the form de = 
—pd(l/D) + Tdb with pressure p(D, b) and temperature T(D, b). 

This action principle fits into the general theory with the electromagnetic field 
variables playing the role of additional configuration variables which are not acted 
on by the particle relabelling group. They obey the usual Euler-Lagrange equa- 
tions, coupled to the Euler-Poincare variables through the Lagrangian. In other 
words, the primitive unreduced Lagrangian in this case is of the abstract form 
L : TG x V* x TQ x TC -> R in which G, the fluid particle relabelling group, 
acts trivially on the Maxwell field variables Q and the gauge field 9 G C. Note that 
the Lagrangian in equation ( [7.26 ) is invariant under translations of 9, as well as 
under the electromagnetic gauge transformations, 



A + VC, $ -»■ $ - dC/dt, 6 



C, 



(7.27) 



for an arbitary function £ of x and t. 

We now take variations of the action. The variation of S re( j in equation (7.26) 
may be written using the definitions of E and B, and the abbreviated notation 
c = A- v — <J> + a), as 



(56 rod 



dtd 3 x 



D (v + cA + cV8) ■ (5v + 



1 



v| 2 + ic 2 -e-^ UD 
1 2 D 



Ot 



dE 
~dt 



+ V • cDv 89 



DT 5b + ( cDv + 



dcD 



VxB -SA + (-cD + V • E) <5$ 



(7.28) 



where terms arising from integration by parts vanish for the natural boundary con- 
ditions given by 



v • n = 0, E • n = 0, and n x B = on the boundary, 



(7.29) 
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and for variations 5g(t) of g(t) vanishing at the endpoints. Stationarity of the action 
S re d in ( |7.26| ) under variation of the gauge field 9 gives the conservation law, 

d \ 3D 

— + v • V j c = via the continuity equation, + V • (Dv) = 0. (7.30) 



Hence, we may set c = q/m in equation (7.28) and then acquire the two Maxwell 
equations with sources from stationarity of the action S re( j under variations of A 
and Once the flow velocity v is known, the relation c = q/m determines the 
gauge function 9 by "quadrature" , from the definitions of c and u> as 

f)9 n 

u= — + v-V0 = -^ + *-A-v. (7.31) 
at m 

The remaining variations of © re d i n { v > D, b} for the Euler-Poincare dynamics 
collect into the Kelvin-Noether form of equation ( |7.9| ) as 

d \ ( 1 61 \ 1 61 „ , 1 61 „, 51 , 

7T + v • V — + — -^-Vv J + — — V6 - V— = . (7.32) 
dt J\D6vJ D5vi D6b 6D v ; 

Specifically, we have 

+ v • V^j (v + cA + cV9) + {vj + cAj + c9 d )Vv j 

-JVb- vQ|v| 2 + ^c 2 -e-|) =0. (7.33) 

Using the fundamental vector identity of fluid dynamics in three dimensions, 

(b • V)a + OjVy = -b x (V x a) + V(a • b) , (7.34) 

with, in this case, b = v and a = D~ l 5l/5v, casts the Euler-Poincare equation 
( 7.32] ) into its equivalent "curl" form, 



d ( 1 51 \ f I 61 \ I 51 / 1 51 51 \ 



Similarly, applying the same vector identity with b = v and a = c(A + V0) in the 
Maxwell fluid motion equation ( |7.33| ) yields, 

dv . t— r\ 1„ / OA ^ . ^ 

+ (v • V)v + — Vp = c — V$ + vx(VxA ; 



dt D V dt 

= — (E + v x B) , (7.36) 

m 

where we have used the thermodynamic first law. Thus we find the Maxwell fluid 
motion law - the first among the equations in ( 7.24j ) - after setting c = q/m and 
using the definitions of the electromagnetic fields E and B in terms of the potentials 
A and <£. 
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Theorem 7.1 (Kelvin circulation theorem for the Maxwell fluid.) By the MF 

motion equation ( \7. 3(\ ) and the thermodynamic first law, we have 



dJ_ 
~dl 



Tdb , (7.37) 



where the circulation integral is given by 

1= I (v + — A) • dx, 

H m 

for a curve jt which moves with the fluid velocity v. 



Proof. The proof is the same as for Theorem p.2[ although it is also immediately 
seen from the motion equation ( 7.33 ) after substituting D _1 5//5v = v + c(A + W) 
and D^Sl/Sb = -T. ■ 

Corollary 7.2 (Potential vorticity convection for the Maxwell fluid.) Stokes' 
theorem, advection of specific entropy b and the continuity equation together imply 
convection of potential vorticity for the adiabatic Maxwell fluid, 

^ + v-Vg = ^ = with q = — V6-curl(v + ^-A). (7.38) 
at dt Dm 



Remark. The equation dq/dt = for convection of potential vorticity for a general 
Lagrangian with dependence l(y,b,D), with 

5 = _L Vi ,.curl(i£), (7.39) 

may also be proven directly from the "curl" form of the Kelvin-Noether equation 
( 7.35| ) in three dimensions, by taking the scalar product of its curl with V6 and 
applying the continuity equation for D. 



Alternative interpretations of the Maxwell fluid formulation. Note that 
the first line of the Euler-Poincare motion equation ( 7.36| ) for Maxwell fluids persists, 
when the electromagnetic energy terms are dropped from the MF Lagrangian in 
equation ( [7.26| ), to give the action 



6 



red 



dt I 



dt d 6 x 



-D|v| 2 + -Z?(A-v 
2 1 1 2 V 



$ + cD) 2 - De(D, b)) . (7.40) 



The Euler-Poincare equation which results from this action is 



dv . 1„ 
_ + (v.V)v + -V P 



dA 

-— - v$ + v x (v x a; 



(7.41) 



The Kelvin Theorem 7.1 and its potential vorticity Corollary |7.2| also persist for the 
dynamics derived from this truncated action. The "Lorentz force" terms in equation 
(7.41) in terms of A and $ arise purely from the Kaluza-Klein coupling term — the 



48 



second term in the integrand of the action S re d in equation ( |7.40| ). These "Lorentz 
forces" may be interpreted physically as noninertial forces resulting from having 
moved into a frame of reference with a prescribed velocity given by A(x, t). The 
velocity v then represents fluid flow relative to this noninertial frame. This situation 
reduces to Faraday driving of the fluid (Faraday [1831]), when A(x, t) corresponds 
to a rigid motion of the fluid container. For a simple example, set c = /o; = 0; 
and A = ^z x x. Then cv x (V x A) = v x /qz gives the Coriolis force and v 
corresponds to fluid velocity in a uniformly rotating reference frame with constant 
angular velocity /o . This is the typical situation in geophysical fluid dynamics. 
Alternatively, the right hand side of equation ( |7.41 ) may be interpreted as a 



vortex force arising from a given wave field at the surface of an incompressible fluid 
(for D = 1), as in Craik and Leibovich [1976] (see also Holm [1996b]). The Craik- 
Leibovich equations are formally identical to equation ( [7.41| ) when A is identified 
as the prescribed mean Stokes drift velocity due to the presence of the wave field. 
Fluid motion equations of the same form as ( 7.4l| ) also appear in the generalized 



Lagrangian-mean (GLM) formulation of wave, mean-flow interaction theories (see 
Andrews and Mclntyre [1978a,b]), in which Co is the Doppler-shifted frequency of 
a wave packet interacting with a Lagrangian-mean flow of velocity v, and A is 
the prescribed pseudomomentum per unit mass of the wave. For a discussion of 
self- consistent Lie-Poisson Hamiltonian theories of wave, mean-flow interaction in a 
similar form, see Gjaja and Holm [1996]. 

Geodesic motion and the Kaluza-Klein construction for incompressible 
fluids. Hamilton's principle for the action in the "minimal coupling" form 

6'= fdtl= [ dtd 3 x (-D\v\ 2 + —DA-v-—D<S>-De(D,b)) , (7.42) 
J J \2 m m J 

yields the same Euler-Poincare equation ( |7.41|) as results from Hamilton's principle 
for the Kaluza-Klein action in equation ( |7.40| ). Thus, we see that by introducing the 
auxiliary gauge field, 6, the Kaluza-Klein construction renders the minimal coupling 
form of the action for fluid dynamics quadratic in the velocity, while preserving its 
corresponding Euler-Poincare equation. The Kaluza-Klein construction for charged 
particle mechanics is discussed in Marsden and Ratiu [1994]. The historical refer- 
ences are Kaluza [1921], Klein [1926] and Thirry [1948]. 

In the incompressible case, when the Kaluza-Klein action is taken to be 

S red = j dtl = J dtd 3 x Qd| v| 2 + ^D(A • v - $ + Co) 2 - p(D - 1)^ , (7.43) 

for arbitrary prescribed functions A and the resulting Euler-Poincare equation 
(i.e., equation ( [7.41 ) with D = 1) represents geodesic motion on the group of volume- 



preserving difleomorphisms with respect to the conserved kinetic-energy metric given 
by 

||v|| 2 = / d 3 x Q|v| 2 + i(A-v-$ + cD) 2 j , (7.44) 
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where c = A • v — + u; is an advected quantity, dc/dt = 0. This observation ex- 
tends the geodesic property of incompressible ideal fluid flows established in Arnold 
[1966a] to the case of incompressible Maxwell fluid flows, as well as to the case of 
incompressible ideal fluid flows in an arbitrarily moving reference frame. From the 
Euler-Poincare point of view, this extension enlarges the particle relabelling group 
G from the group of diffeomorphisms to the group of automorphisms of the single 
particle Kaluza-Klein bundle. The total system for the incompressible Maxwell fluid 
flows is then geodesic motion on the product of this automorphism group with the 
Maxwell fields themselves. We believe that a similar extension may be involved in 
the results of Ono [1995a, 1995b]. 

We should also remark that when one has equations in geodesic form, one can 
make use of all the attendant geometry to obtain additional interesting results. 
Examples of this applied to questions of stability and conjugate points are given in 
the works of Misiolek listed in the references. 

8 Approximate Model Fluid Equations which Preserve 
the Euler-Poincare Structure 

The preceding section demonstrates the applicability of the Euler-Poincare theorem 
for ideal continua when the equations of motion are given. Here we discuss approx- 
imate fluid models which preserve the Euler-Poincare structure, and are obtained 
by making asymptotic expansions and other approximations in Hamilton's principle 
for a given set of model equations. As examples, in this section we first discuss the 
derivation of the quasigeostrophic approximation in geophysical fluid dynamics from 
an approximation of Hamilton's principle for the rotating shallow water equations. 
Next, we discuss the Boussinesq approximation for dispersive water waves in one 
dimension. As an example of the type of "bonus" which may appear in making 
simplifying approximations while preserving mathematical structure, we derive the 
integrable Camassa-Holm equation (Camassa and Holm [1993], Camassa, Holm and 
Hyman [1994], Alber et al. [1994, 1995, 1997]), by making asymptotic approxima- 
tions in the Hamilton's principle for the Boussinesq equations. The Camassa-Holm 
equation in one dimension is a completely integrable partial differential equation for 
dispersive water waves that was actually discovered by making structure preserving 
approximations of this type. This equation turns out to describe geodesic motion 
on the group of diffeomorphisms of either the real line or the periodic interval, with 
metric given by the H 1 norm of the velocity. We also derive a multidimensional ana- 
logue of the one-dimensional Camassa-Holm equation by invoking the n-dimensional 
version of this geodesic property. There are also other potential advantages of mak- 
ing structure preserving approximations, e.g., for numerical integrations. However, 
discussion of these other advantages is deferred to another place. (See Marsden 
and Wendlandt [1997], Wendlandt and Marsden [1997] and Marsden, Patrick and 
Shkoller [1997] for recent advances in this direction.) 
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Rotating shallow water dynamics as Euler Poincare equations. We first 
consider dynamics of rotating shallow water (RSW) in a two dimensional domain 
with horizontal coordinates x = (2:1,2:2). RSW motion is governed by the following 
nondimensional equations for horizontal fluid velocity v = (^1,^2) an d depth D, 

d d D 

e-v + /^xv + V^ = 0, — + V-£>v = 0, (8.1) 
at at 



with notation 



d fd „\ , , f D — B 
— = I — + v • V I and ip 



dt \dt / \ 



These equations include variable Coriolis parameter / = /(x) and bottom topogra- 
phy B = B(x). 



The dimensionless scale factors appearing in the RSW equations (8.1) and ( |8. 2| ) 
are the Rossby number e and the rotational Froude number J 7 , given in terms of 
typical dimensional scales by 

e = -^ r <l and T = ^! = 0(1) . (8.3) 
JqL gB 

The dimensional scales (Bq, L, Vo, /o, <?) denote equilibrium fluid depth, horizontal 
length scale, horizontal fluid velocity, reference Coriolis parameter, and gravitational 
acceleration, respectively. Dimensionless quantities in equations Q8.1| ) are unadorned 
and are related to their dimensional counterparts (primed), according to 

v' = V v, x' = Lx, f = t, f = M, 



\V , 

B' = B B, D' = B D, and D' - B' = B (D - B). (8.4) 

Here, dimensional quantities are: v', the horizontal fluid velocity; D', the fluid 
depth; B', the equilibrium depth; and D' — B' , the free surface elevation. 

For barotropic horizontal motions at length scales L in the ocean, say, for which 
T is order 0(1) - as we shall assume - the Rossby number e is typically quite small 
(e ^ 1) as indicated in equation (|3.3D . Thus, e < 1 is a natural parameter for 
making asymptotic expansions. For example, we shall assume |V/| = 0(e) and 
|V-B| = O(e), so we may write / = 1 + e/i(x) and B = 1 + e-Bi(x). In this scaling, 
the leading order implications of equation ( |8.1| ) are v = z x V*0 and V • v = 0. This 
is geostrophic balance. 

A simple calculation using equation ( 7.1 2| ) shows that the RSW equations ( |3,lD 



arise as Euler-Poincare equations from Hamilton's principle with action 6rsw> 



Srsw = / di/R,sw = dt dxidx 2 



Dv . R(x) _Pz|)! + ! fl | v p 



(8.5) 



2eT 2 

where curlR(x) = /(x)z yields the prescribed Coriolis parameter. The RSW equa- 
tions (^^) themselves can be regarded as being derived from asymptotics in Hamil- 
ton's principle for three dimensional incompressible fluid motion, see Holm [1996a]. 
However, this viewpoint is not pursued further here, as we proceed to describe the 
relation of RSW to the quasigeostrophic approximation of geophysical fluid dynam- 
ics. 
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Quasigeostrophy. The quasigeostrophic (QG) approximation is a useful model 
in the analysis of geophysical and astrophysical fluid dynamics, see, e.g., Pedlosky 
[1987]. Physically QG theory applies when the motion is nearly in geostrophic 
balance, i.e., when pressure gradients nearly balance the Coriolis force in a rotating 
frame of reference, as occurs in meso- and large-scale oceanic and atmospheric flows 
on Earth. Mathematically the simplest case is for a constant density fluid in a 
planar domain with Euclidean coordinates x = (aTi,^)- QG dynamics for this case 
is expressed by the following nondimensional evolution equation for the stream- 
function tp of the incompressible geostrophic fluid velocity v = z x Vip, 

Here A is the Laplacian operator in the plane, J- denotes rotational Froude number, 
[a, b] = d(a,b) /d(xi, X2) is the Jacobi bracket (Jacobian) for functions a and b 
defined on M? and j3 is the gradient of the Coriolis parameter, /, taken as / = l+f3x2 
in the /?-plane approximation, with constant 0. (Neglecting (3 gives the /-plane 
approximation of QG dynamics.) The QG equation ( |8.6| ) may be derived from an 
asymptotic expansion of the RSW equations ( |8.1|) by truncating at first order in the 
Rossby number, cf. Pedlosky [1987]. Equation Q8,6|) may be written equivalently in 



terms of the potential vorticity, q, as in equation (7.3- 



|| + v • Vg = 0, where q = Aijj - T$ + / for QG. (8.7) 

This form of QG dynamics expresses its basic property of potential vorticity con- 
servation on geostrophic fluid parcels. 

The QG approximation to the RSW equations introduces "quasigeostrophic par- 
ticles" which move with geostrophic velocity v = z x X7ip and, thus, trace the 
geostrophic component of the RSW fluid flow. These QG fluid trajectories are de- 
scribed as functions of Lagrangian mass coordinates i = (^1,^2) given by x(€, t) in 
the domain of flow. 



Hamilton's principle derivation of QG as Euler Poincare equations. As 

in Holm and Zeitlin [1997], we consider the following action for QG written in the 
Eulerian velocity representation with the integral operator (1 — JFA -1 ), 



6 



red 



dtl 



dt I dx\dx2 



Dv(l- .FA _1 )v + Dw ■ R(x) - ip(D - 1) 



This choice can be found as an asymptotic approximation of the RSW action Grsw 
in equation (8.5), in the limit of small wave amplitudes of order 0(e 2 ) and constant 



mean depth to the same order, when the wave elevation is determined from the 
fluid velocity by inverting the geostrophic relation, v = z x Vf/>. The variational 
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R + e 

e 



T A~\ - Z-A-^Dv) 



2 2L> 
v(l-.FA _1 )v + v ' R - 



(D-l) 



derivatives of the reduced Lagrangian 6 rc d at fixed x and t are 

Id 2 

where we have used the symmetry of the Laplacian operator and assumed no con- 
tribution arises from the boundary when integrating by parts. For example, we 
may take the domain to be periodic. Hence, the Euler-Poincare equation ( 7.12| ) for 
action principles of this type and the fundamental vector identity ( 7.34| ) combine to 
give the Eulerian QG "motion equation" , 

d 

:(1 - ^'A" 1 )v - v x curl (e(l - TA^v + R) 



upon substituting the constraint D 
equation yields 

dq 



0, (8.10) 
1, imposed by varying ip. The curl of this 



— + v • Vq + gV • v 

ot 



0, 



where the potential vorticity q is given by 

q = ez • curl (1 - ^*A -1 )v + / = e(A^ - Tip) + / . 

with 

/ = z • curlR = 1 + (3x2, 



ill 



U2) 



5.13) 



and (3 is assumed to be of order 0{e). The constraint D = 1 implies V • v = (from 
the kinematic relation dD/dt + V ■ Dv = 0) and when v = z x Vip is substituted, the 
equation for q = Aip — Tip + / yields the QG potential vorticity convection equation 



(8.7). Thus, the QG approximation follows as the Euler-Poincare equation for an 
asymptotic expansion of the action for the RSW equations when the potential energy 
is modelled by inverting the geostrophic relation. The same QG equation follows 
upon recasting the action fl3.8| ) in the Kaluza-Klein form ( |7.43j ) for incompressible 
fluids, 



6 



red 



dtl 



dt dx\dx2 



-Dv • (1 

2 v 



^A _1 )v + i-D(R • v + Cjf - ip(D - 1 



1 



~.\2 



J.14) 



where Q is defined as u> = d9/dt = 89/ dt + v • V6 for the gauge field 6, as in the 
case of the Maxwell fluid. Thus, the QG motion equation (|8.10| ) with the beta- 
effect (included in R) describes geodesic motion on the group of area-preserving 
diffeomorphisms with respect to the conserved kinetic-energy metric given by 



dx\dx2 



|v(l-^A- 1 )v+i(R.v + 1 i) 2 



5.15) 
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where d = R • v + ui is an advected quantity, dc' /dt = 0. This observation from the 
Euler-Poincare viewpoint confirms the geodesic interpretation of the QG equations 
for motion in the /3-plane established in Zeitlin and Pasmanter [1994]. 



ID Boussinesq dispersive shallow water equations. For one dimensional 
shallow water motion with prescribed mean depth B{x) we choose the following 
action 



S re d 



/*!-/ 



dt dx 



l -Dv 2 + ^(Dv) 



9 



(D-B(x)) 2 



U6) 



in which g and a 2 are constants and subscript x denotes partial derivative. The 
second term, proportional to a 2 , represents the kinetic energy due to vertical motion. 
The last term is the potential energy. Recall that the surface elevation h = D — B{x) 
satisfies dh/dt = —{Dv) x for shallow water dynamics in one dimension. Thus, the 
last two terms are analogous to the Lagrangian 



a 2 / 



dhy 
2 K~dt ) 



9 -h 2 
2 



This is the Lagrangian for a harmonic oscillator whose displacement is given by the 
surface elevation h and whose natural frequency is y/g/a. For the choice of action 
in equation ( 8.1q* ) and for boundary conditions such that v x — ► as |x| — ► oo, the 
variation is 



S& 



red 



dt dx 



D(v-a 2 (Dv) xx )5v + 



a'v{Dv) xx -gh)8D 



3.17) 



The Kelvin-Noether form of the Euler-Poincare equations in ( |7 . 1 2| ) then gives 
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dt + tv )\D5v dX 



<w) 



d 

— (v - a 2 (Dv) 3 



+ vv x + gh a 



dx. (8.18) 



Upon inserting the one-dimensional continuity equation into the a 2 term and rear- 
ranging slightly, we find the system of equations 



dv , id 2 h x 

- Ft +vv x +gh x + a-^ r 



dh 

0, — + (hv + B(x)v) x 



0. 



3.19) 



the second of which is just a restatement of continuity. According to Whitham 
[1974] these equations were favored by Boussinesq, who first formulated them by 
using the method of asymptotic expansions. Here we see that the Boussinesq shallow 
water equations ( 8.19 ) are also Euler-Poincare equations on Diff(M) derived from 



the action ( |8.16| ). The term proportional to a 2 in equation ( 8.19] ) arises from the 
kinetic energy due to vertical motion in the action (8.16) and produces the wave 



dispersion responsible for solitary wave solutions of these equations. 
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ID Camassa Holm equation for peakons. In the limit that the potential 
energy gh 2 jl is negligible compared to the kinetic energy (e.g., for weak gravity, or 
small surface elevation), we may ignore the last term in the action ( 8.16| ) for the 



Boussinesq equations in one dimension, set D = B(x) in the other terms and rescale 
to c? = 1, conforming to the notation in Camassa and Holm [1993]. We thereby 
obtain the following simplified expression for the shallow water action in this regime, 

6 red = J dtdxB(x) fi v 2 + ^A _ (g 20) 

For this action, when we also assume B{x) = 1 (for constant bottom topography) 
Hamilton's principle implies simply 

= J6 rcd = jdtl = J dtdx (v - v xx )5v , (8.21) 

for vanishing boundary conditions for v x on the real line as \x\ — > oo. Hence, 
51 /5v = v — v xx and the basic Euler-Poincare equations for this case reduce to 

This is the k = case of the completely integrable partial differential equation 
derived by Camassa and Holm [1993], 

d 

777 (v - v xx ) + 2kv x = -3vv x + 2v x v xx + vv xxx . (8.23) 
at 

For k = 0, this equation admits 'peakon' solutions. The peakons are solitons which 
interact elastically and possess a peak, at which the derivative v x reverses sign. The 
simplest case is the single peakon, which is a solitary travelling wave solution given 
by v(x,t) = Co exp — \x — cot\, with a constant wave speed cq. The multi-peakon 
solutions of the Camassa-Holm equation are obtainable from its associated Lax pair 
and linear isospectral problem, as shown in Camassa and Holm [1993]. 

Being basic Euler-Poincare, equation ( 8.22; ) describes geodesic motion. Camassa 



and Holm [1993] note that the integrable dynamics of N peakons interacting nonlin- 
early via equation ( |8.22j ) reduces to finite dimensional geodesic motion on a manifold 
with N corners. This geodesic property persists to infinite dimensions and although 
the equation was originally intended to be an approximation of shallow water mo- 
tion, it turns out equation ( 8.22|) is also the geodesic spray equation for motion on 



the group of diffeomorphisms of the real line with metric given by the H 1 norm of 



v, see Kouranbaeva [1997]. The k / case of the Camassa-Holm equation ( 8.23 ) 
may be obtained formally by shifting (v — v xx ) by k in equation ( |8.22j ) and retaining 
homogeneous boundary conditions for (v — v xx ) as |x| — > oo. The corresponding 
statement about geodesic motion for however, is rather more technical than 

for k = and involves the Gel'fand-Fuchs co-cycle and the Bott-Virasoro group, see 
Misiolek [1997] for details. See Alber et al. [1994, 1995, 1997] for discussions of the 
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periodic solutions of the Camassa-Holm equation and a related integrable shallow 
water equation in the Dym hierarchy, 



2kv x 



a 

a? 



+ 2v x v xx + vv 

XXX ■ 



(8.24) 



This equation is the "high wave number limit" of the Camassa-Holm equation ( 8. 23] ). 



Higher dimensional Camassa-Holm equation. As we have seen, the Camassa- 
Holm (CH) equation in one dimension describes geodesic motion on the diffeomor- 
phism group with respect to the metric given by the H 1 norm of the Eulerian 
fluid velocity. Thus, a candidate for its n-dimensional incompressible generalization 
should be the Euler-Poincare equation that follows from the Lagrangian given by 
the H 1 norm of the fluid velocity in n dimensions, subject to volume preservation 
(forn^l), 



6 



red 



dt I 



D 



dt I d n x —{viv 1 + aX»') -p(D - 1) 
M 2 



5.25) 



where we have restored the length-scale, or aspect-ratio parameter, a. Varying this 
action at fixed x and t gives 



5e 



red 



dt I d n x 

M 

+ a : 



v % + a 2 vfv l j)5D - (D — l)Sp + [Dvi - a 2 dj{Dvf)) 8v 



I dtl d n ~ X : 
J JdM 



x hj(Dvf 5v l ) , 



whose natural boundary conditions on dAi are 

v • n = and (n • V)v || n, 



(8.26) 



5.27) 



where || denotes "parallel to" in the second boundary condition, which of course 
is not imposed when a 2 is absent. (Recall that <5v in equation (|S.26) is arbitrary 



except for being tangent on the boundary. This tangency, along with the second 
condition in equation ( |S.27 ) is sufficient for the boundary integral in equation (8.26) 
to vanish.) By equation (7J3) or ( 7.32 ), the Euler-Poincare equation for the action 



S re d in equation ( 8.25 ) is 
d 



dt 



+ v • V (v - a 2 Av) + (vj - a 2 Avj)Vv J - V 



+ 



a 



|Vv| 



P 



0. 

3.28) 



where (Vv) 



dv l /dxi, |Vv| 



tr(Vv • Vv T ) and superscript 



(•) T denotes transpose. We have also used the constraint D = 1, which as before 
implies incompressibility via the continuity equation for D. Requiring the motion 
equation ((3.28D to preserve divv = implies a Poisson equation for the pressure 
p with a Neumann boundary condition, which is obtained just as in the case of 
incompressible ideal fluid dynamics by taking the normal component of the motion 
equation evaluated at the boundary. 
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Properties of the Camassa Holm equation. Since the CH action 6 re d in. 
^25| ) is translation invariant, the Noether theorem ensures the CH equation ( 8.28 ) 



conserves a momentum. In fact, by the stress tensor formulae Q7.6| )-( |7T8| ), equation 

T( . (8.29) 



).2q ) may be rewritten as 

drrij d 



dt dxi 

In this case, the momentum density m;, z = 1,2,3 defined in equation ( |7.7[ ) is given 
by 



1 5v l 



= Vi-a 2 Avi, (8.30) 

D=l 



and the stress tensor T- defined in equation (fflq ) is given by 

Tl = ( Vi - a 2 A Vi y - a 2 ufyjjf + tip . (8.31) 



Thus, equation ( 8.29j ) implies conservation of the total momentum, M = J M md 



x. 



provided the normal component of the stress tensor T- vanishes on the boundary. 

Since the CH equation (|8.28j ) is Euler-Poincare, it also has a corresponding 
Kelvin- Noether circulation theorem. Namely, cf. equation ( p. 32 ), 



d_ 

It 



(v - a 2 Av) • dx = 0, (8.32) 



for any closed curve jt that moves with the fluid velocity v. This expression for the 
Kelvin-Noether property of the CH equation in 3D is reminiscent of corresponding 
expressions in wave, mean-flow interaction theory. This correspondence suggests a 
physical interpretation of the a 2 term in the Kelvin-Noether circulation integral as 
a Lagrangian mean closure relation for the pseudomomentum of the high frequency 
(i.e., rapidly fluctuating, turbulent) components of the flow. In this interpretation, 
a corresponds to the typical length scale at which these high frequency components 
become important. See Holm, Foias and Titi [1998] for more discussion of using the 
3D CH equation (|8.28j) as the basis for a turbulence closure model. 



In three dimensions, we may use the vector identity ( 7.34 ) to re-express the CH 



motion equation ( ^.28[) in its "curl" form, as 
d 



dt 



1 - a 2 A)v - v x (V x (1 - a 2 A)v) 



+ V ( v (l-a 2 A)v-i|v| 2 -y|Vv| 2 +p j - (). (*.:«) 



The inner product of v with this equation then implies conservation of energy, 

E= \ I d 3 x (v • (1 - a 2 A)v) = \ [ d 3 x (|v| 2 + q 2 |Vv| 2 ) , (8.34) 
2 Jm ^ J 



upon integrating by parts and using the boundary conditions ( 8.27|) . Naturally, this 



energy is also conserved in n dimensions. In fact, Legendre transforming the action 
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(|8.25| ) gives the following Hamiltonian (still expressed in terms of the velocity, 
instead of the momentum density m = Sl/5v), 



H 



[ d n x (|v| 2 + a 2 |Vv| 2 ) +p(D- 1)1 . (8.35) 
Jm l 2 J 



Thus, when evaluated on the constraint manifold D = 1, the Lagrangian and the 
Hamiltonian for the CH equation coincide in n dimensions. (This, of course, is not 
unexpected for a stationary principle giving rise to geodesic motion.) 



The curl of the 3D Camassa-Holm motion equation ( 8.33 ) yields 



d 

— q = q • Vv — v • Vq = [v, q], where q = curl(v — a 2 Av) , (8.36) 

and we have used incompressibility and commutativity of the divergence and Lapla- 
cian operators. Thus, v is the transport velocity for the generalized vorticity q and 
the "vortex stretching" term q- Vv involves Vv, whose I? norm is controlled by the 



conservation of energy in equation (8.34). Boundedness of this norm will be useful 
in future analytical studies of the 3D Camassa-Holm equation; for example, in the 
investigation of the Liapunov stability properties of its equilibrium solutions. 

3D periodic CH motion. In a three dimensional periodic domain, the conserved 
energy E in equation ( |3.34| ) may also be expressed as 

E = - [ d 3 x (|v| 2 + a 2 |curlv| 2 ) , (8.37) 

upon integrating by parts and using divv = 0. Thus, in the 3D periodic case, the 
CH energy E may be interpreted as the sum of the kinetic energy and the enstrophy 
(i.e., the 1? norm of vorticity) of the Euler fluid. 



The inner product of the generalized vorticity q with the motion equation (8.33) 
implies conservation of helicity, for three dimensional periodic motion. Namely, 
the quantity 



A 



/ d 3 x (1 - a 2 A)v • curl(l - a 2 A)v (helicity), (8.38) 
Jm 



I M 

is also a constant of motion for three dimensional periodic CH motion. 



^From the CH vorticity equation ( |8.36j ), we see that steady 3D solutions of 
the CH equation (denoted with subscript e for "equilibrium") are characterized by 
the vector-field commutation relation [v e ,q e ] = 0. Thus, the velocity of a steady 
CH flow v e generates a volume preserving diffeomorphism that leaves invariant its 
corresponding steady generalized vorticity q e . For example, the CH Beltrami 
flows for equation ( 8,33j ) are characterized by v e = /iq e , for a constant [i. The 



CH Beltrami flows verify the invariance property; hence, they are steady. These 
steady solutions are also critical points of the sum E + A/2/i of the energy E in 
equation ( |3.37 ) and 1/2^ times the conserved helicity in equation ( |8.38j ). Hence, 



they are relative equilibrium solutions of the CH equation. The CH Beltrami flows 
are divergenceless vector eigenfunctions of the product of the curl operator and the 
Helmholtz operator, (1 — a 2 A). They are the CH analogues of "A-B-C flows" for 
the ideal Euler fluid. 
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Constitutive properties of the CH "fluid." Physically, conservation of the 
energy in equation ( |3.37| ) means that the CH fluid can exchange energy between its 
translational, and its rotational and shear motions. One may ask, what constitutive 
relation describes such a fluid? 

One may verify directly that the 3D Camassa-Holm equation ( |8.28| ) in Cartesian 
coordinates implies the following formula for the geodesic spray form of the CH 
equations in 3D: 



(l-n-'A)(^+vV) r, ■= rr'lA,;,)!,;;-.-') 



d 



a 



P ~ Y V * V '* ) S i + 2a2y3 ,k v f 



5.39) 



In vector notation, this is 



( L-n 2 A)(^ + vV) v 



a 2 (Av) x curlv 

,2 



div 



P ~ y l Vv l ) 1 + 2a z Vv • Vv J 



5.40) 



where I is the unit tensor. Viewed this way, the CH fluid acceleration dv/dt is non- 
local, nonlinear, and (as we know from the general Euler-Poincare theory) geodesic. 
Holm, Foias and Titi [1998] show, among other things, that the geodesic spray 
equation ( |3.40D may be rearranged into a viscoelastic constitutive relation of Jef- 
freys type. 



Discussion. The essential idea of the CH equation is that its specific momentum 
(i.e., its momentum per unit mass) is transported by a velocity which is smoothed 
by inverting the elliptic Helmholtz operator (1 — a 2 A), where a corresponds to the 
length scale at which this smoothing becomes important, i.e., when it becomes of 
order O(l). When the smoothing operator (1 — a 2 A)" 1 is applied to the transport 
velocity in Euler's equation to produce the CH equation, its effect on length scales 
smaller than a is that steep gradients of the specific momentum function tend not 
to steepen much further, and thin vortex tubes tend not to get much thinner as they 
are transported. And, its effect on length scales that are considerably larger than a 
is negligible. Hence, the transport of vorticity in the CH equation is intermediate 
between that for the Euler equations in 2D and 3D. As for Euler vorticity, the curl 
of the CH specific momentum is convected as an active two form, but its transport 
velocity is the smoothed, or filtered CH specific momentum. 

The effects of this smoothing or filtering of the transport velocity in the CH 
equation can be seen quite clearly from its Fourier spectral representation in the 
periodic case. In this case, we define nik as the k-th Fourier mode of the specific 
momentum m = (1 — a 2 A)v for the CH equation; so that nik = (1 + a 2 |k| 2 )vk. 
Then the Fourier spectral representation of the CH equation for a periodic three- 
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dimensional domain is expressed as 




1 + a 2 |p| 2 

p+q=k ^ IHI 



(q x m q ) = 0, (8.41) 



where IT is the Leray projection onto Fourier modes transverse to k. (As usual, the 
Leray projection ensures incompressibility.) In this Fourier spectral representation 
of the CH equation, one sees that the coupling to high modes is suppressed by the 
denominator when l+a 2 |p| 2 3> 1. Consequently, when |p| > 0(l/a), the smoothing 
of the transport velocity suppresses the development of higher modes |k|>0(2|p|). 
And, it also suppresses the "stochastic backscatter" from higher modes to lower 
ones, |k| = 0(1). Thus, thinking of "interaction triangles" among the modes, one 
sees that both the very acute and very obtuse p + q = k triangles are suppressed, 
leaving active only the nearly equilateral ones, when |p| > 0(l/a). Hence, the 
CH smoothing of the transport velocity suppresses both the forward and backward 
cascades for wave numbers |p| > 0(1/ a), but leaves the Euler dynamics essentially 
unchanged for smaller wave numbers. As we have seen, the result is that the vortex 
stretching term in the dynamics of q = curl m is mollified in the CH model and so 
the vortices at high wave numbers will tend to be "shorter and fatter" than in the 
corresponding Euler case. 

When the kinetic energy terms are neglected relative to the gradient velocity 
terms, the CH action S re d in ( |B.25 ) becomes 



©red = f dtl= I dt [ d n x 
J J J M 



jv^vf -p(D-l) 



(8.42) 



whose Euler-Poincare equation in 3D implies the following, 

^Av - v x (V x Av) + V ^v • Av + ^|Vv| 2 - p^j = . (8.43) 

This equation is the "high wave number limit" of the 3D Camassa-Holm equation 
33| ). Scale invariance is restored in this limiting equation and its corresponding 



group invariant (e.g., self-similar) solutions may be illuminating. 



The Riemannian CH equations. One can formulate the CH equations on a 
general Riemannian manifold, possibly with boundary. Although this will be the 
subject of future papers, we will make some comments about some of the features 
(some of them conjectural) here. 

We start with a smooth, oriented, compact Riemannian manifold M, possibly 
with a smooth boundary. We first define the group DiffcH to be the group of 
diffeomorphisms r] : M — > M of class H s , where s > (ra/2) + 2 with the boundary 
condition that the tangent map Tn : TM — > TM takes the outward normal direction 
to the boundary dM at a point x £ M to the outward normal direction at the point 
r](x). We first conjecture that this group is a smooth manifold and is a Lie group 
(in the same sense as in Ebin and Marsden [1970]) with Lie algebra the set of vector 



60 



fields v on M which are tangent to the boundary of M and that satisfy the boundary 
condition 



(V n v,u) = S(u,v) 



for all vectors u tangent to the boundary. Here S(u, v) is the second fundamental 
form of the boundary. This condition on the boundary is the CH analogue to 
the condition of parallel to the boundary in the case of the Euler equations. The 
condition comes about by differentiation of the condition on rj in the definition of 
DiffcH using a routine calculation. 

Notice that the boundary conditions are different from those previously (see 
equation ( |8.27j ). This appears, however, to be needed for the group theoretic version 
of the equations. Now we put a right invariant Lagrangian on DiffcH which, at the 
identity, is given by 



where Vt> is the covariant derivative of v, where dfi is the Riemannian volume 
element and where dA is the area element of the boundary. We are using terminology 
appropriate to the case in which M is three dimensional, but of course there is no 
restriction on the dimension of M. Also, ||Vt> || 2 denotes the norm in the sense of 
the full Riemannian contraction of the tensor Vv. The associated Laplace operator 
is usually called the rough Laplacian. 

At this point there are some choices one can make. One can use a different H 1 
metric built out of thinking of v as a one form and using the d and 5 operators 
and the corresponding Laplace deRham operator. This leads to a slightly different 
system in general, but one that has similar analytical properties. 

The Riemannian CH equations are, by definition, the Euler-Poincare equa- 
tions for this group and this Lagrangian. The boundary term in the Lagrangian is 
designed to make the boundary conditions in the resulting equations come out agree- 
ing with those for the Lie algebra of the group DiffcH- Apart from the boundary 
conditions, the resulting equations agree with the ones we developed in Euclidean 
space, but in general one replaces the Laplacian with the rough Laplacian. Note 
that since the Lagrangian is quadratic in v, the equations on DiffcH axe geodesic 
equations (possibly with respect to an indefinite metric). 

Conjecture 8.1 As in the case of the Euler equations (Ebin and Marsden [1970]), 
the geodesic spray of the Riemannian CH equations is smooth if s > (n/2) + 2. 

This conjecture is based on a direct examination of the expression for the spray of 
the Riemannian CH equations. If this is true, then other analytic things, including 
results on the limit of zero viscosity (or viscoelasticity) also hold. We also note that 
because energy conservation involves a stronger norm than in the Euler equations 
for ideal flow, one expects other analytic properties of the Riemannian CH equation 
to be improved. This would include results on stability and long time existence. 



L(v) = \ [ (IM| 2 + 

1 JM 
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Another consequence of this would be that the spray of the incompressible Rie- 
mannian CH equations would also be smooth. This follows since the projection map 
is smooth and the fact that the spray of the incompressible equations is given by the 
composition of the spray of the compressible ones and the tangent of the projection 
map. (These facts are proved in Ebin and Marsden [1970]). 

As we have mentioned, all of these things will be explored in detail in other 
publications. 



2D Camassa Holm equation. In two dimensions, the curl of the Euler-Poincare 
motion equation ( 8.33| ) produces a scalar relation for potential vorticity convection, 
namely, 

^ + v • Vq = 0, where q=(l- a 2 A)A^> for 2dCH. (8.45) 

In terms of the stream function tp, with v = z x Vip, the boundary conditions (8.27) 
in two dimensions with s = z x ri become 

ip = const and n • VV^ • n = Aip = on the boundary. (8.46) 



Potential vorticity convection ( 8.45 ) for the 2dCH equation, combined with incom- 
pressibility and the first boundary condition in (8.27) imply conservation of the 
following quantity, 

C<j> = j d 2 x <S>(q) , (8.47) 

for any suitably well-behaved function <]?. (This is the Casimir, in the Lie-Poisson 
bracket formulation.) Substituting v = z x V^> and using the divergence theorem 
yields an expression for the kinetic energy Lagrangian for the 2dCH equation in 
terms of the stream function, ip. Namely, 

E = ~ J d 2 x (v • (1 - q 2 A)v) = ^Jd 2 x (Vtfj ■ (1 - q 2 A)V^) (8.48) 

= \ ^ ^ W j> ds ^-(1 - a 2 A)^ -^Jd 2 x (V>(1 - a 2 A)A^) , 

where we sum over the connected components of the boundary and use 



constant on the ith component. Thus, solutions of the CH equation ( 8.28j ) in two 
dimensions optimize the integrated product of the stream function and potential 
vorticity, constrained by the circulation of v — a 2 Av on each connected component 
of the boundary. The Lagrange multiplier for this circulation is the corresponding 
boundary component's (constant) stream function. 



Steady CH solutions in two dimensions. Steady solutions of the 2dCH equa- 
tion ( p. 45 ) with boundary conditions ( |8.46 ) satisfy 



z • We x Vg e = J(q e ,ipe) = 0, for q e = (1 - a 2 A)A^ e 



3.49) 
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Thus, steady CH solutions exist when there is a functional relation between the 
potential vorticity q e and its associated stream function tp e . For example, one could 
have q e = F(vp e ) for function F. In particular, the linear steady flows satisfy 

q e = A(l - a 2 A)ip e = -|k| 2 (l + a 2 |k| 2 )V> e . (8.50) 

These are sines and cosines with wave number k for periodic boundary conditions. 
The corresponding fluid velocity for any of these steady 2dCH solutions is found 
from v e = z x VV'e- 

As another example, the point potential-vortex solution centered at z' = 
x' + iy' in the infinite x, y plane has stream function 

jf>(\z -z'\) = \og{\z -z'\)+ K (\z - z'\/a), (8.51) 

where Kq(\z — z'\/a) is the Bessel function of the second kind. This stream function 
satisfies 

q = A(l - a 2 A)V> = 2tt5{\z - z'\), 

where A is the Laplacian operator in the plane and 6{\z — z'\) is the Dirac delta 
function. (The proof uses Alog(|z-z'|) = 2ir5(\z-z'\) and (l-a 2 A)K (\z-z'\/a) = 
2tt5(\z—z'\).) The circular potential-vortex patch solution, q = A(l — a 2 A)ip = 
C = const for \z — z'\ < a and q = for \z — z'\ > a, has stream function 

V> = c( \ Z ~*\ + 1^ for \z-z'\ <a, 

i) = logflz - z'\) + a 2 K (\z - z'\/a) for \z - z'\ > a, (8.52) 

where the constant C is chosen so the velocity is continuous at \z — z'\ = a. The 
interior of this solution is also a uniformly rotating vortex patch. These special 
potential-vortex solutions illustrate the "screening" in the vortex interaction dy- 
namics for the 2dCH equation introduced by the Helmholtz operator (1 — a 2 A), 
which modifies the momentum density of the 2dCH flow relative to the standard 
incompressible Euler equations in the plane. As we have seen, the corresponding 
screening length a is an additional parameter in the model. 



Quasigeostrophic analogue of CH in two dimensions. We extend the nondi- 
mensional QG action principle in equation QB.Sj ) for QG dynamics in a periodic 
two-dimensional domain to include CH a 2 terms, as follows, 



6 



red 



dtl 



dt I dx\dx2 



L2 



Dv ■ (1 - FA' 1 - a 2 A)v + Dv ■ R — ip(D - 1) 



3.53) 



The corresponding Euler-Poincare equation is to be compared to the QG motion 
equation ( |3.10| ). We find, 

d 

e— (1 - J^A" 1 - a 2 A)v - v x curl (e(l - FA~ l - a 2 A)v + R) 



+V (V> + -v • (1 - FA' 1 - a 2 A)vj = . (8.54) 



63 



The curl of this equation yields 

e^ + v-V(eg + /) = 0, (8.55) 
where the potential vorticity q is now given by 

g = z • curl (1 - JFA -1 - a 2 A)v , (8.56) 

and we choose 

/ = z • curlR = 1 + ef3x 2 . (8.57) 

The constraint D = 1 implies V-v = as usual and when v = z x VV 1 is substituted, 
the equation for q = A(l — a 2 /S)ijj — J-tp yields 

^(A(l - a 2 A)^ - + [J,, A(l - a 2 A)V>] + (5^- = . (8.58) 

Here [a, b] = d(a, b)/d(xi, x 2 ) = J (a, b) is the Jacobi bracket (Jacobian) for functions 

a and b defined on the two dimensional domain. Steady solutions of the QG-CH 
equation ( |3.58| ) satisfy 

J(q e + Px 2 ,4> e ) =0, for q e = A(l - a 2 A)^ e - T^ e . (8.59) 



The dispersion relation for plane- wave solutions of equation (8.58) with frequency 
u and wavenumber k is 

Such plane- wave solutions are analogous to Rossby waves in QG. As with QG Rossby 
waves, these plane- wave solutions are nonlinear solutions of equation ( (B.58p in a two 
dimensional periodic domain. If we choose T = 0(1) and a 2 = o(l), then the effect 
of the a 2 term in this dispersion relation is to significantly reduce the oscillation 
frequency and propagation speeds of those waves at wave numbers greater than 
about a -1 . Thus, such short waves are suppressed, and the emerging dynamics of 
this modified QG theory will tend to possess significant activity only at length scales 
larger than a. Apparently the dispersion relation for the dynamics at these larger 
length scales will faithfully approximate the corresponding QG dynamics at these 
scales. This scenario is, of course, consistent with our earlier discussion of the wave 
number dynamics of the CH solutions in three dimensions. 



Appendix 

Left Representation and Right Invariant Lagrangian. There is a version of 
this theorem for right invariant Lagrangians, but with the representation of G on V 
still on the left. The proof is, of course, identical so we shall only state this theorem. 
The set-up is the following: 
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• There is a left representation of Lie group G on the vector space V and G acts 
in the natural way on the right on TG x V*: (v g , a)h = (v g h, h~ l a). 

• Assume that the function L : TG x V* — > R is right G-invariant. 

• In particular, if ao € V*, define the Lagrangian L ao : TG — > R by L ao (v g ) = 
L(v g , ao). Then L ao is right invariant under the lift to TG of the right action 
of G ao on G. 

• Right G-invariance of L permits us to define I : q x V* — > R by 



Conversely, this relation defines for any I : g x V* — > R a right G-invariant 
function L : TG xP-tl. 

• For a curve G G, let ^(t) := ^(f)^(t) -1 and define the curve a(t) as 
the unique solution of the linear differential equation with time dependent 
coefficients d(t) = £(t)a(t) with initial condition a(0) = ao. The solution can 
be equivalently written as a(t) = g(t)ao- 

Theorem 8.2 The following are equivalent: 
i Hamilton's variational principle 



l(v g g ,ga ) = L(v g ,a ). 




(8.61) 



holds, for variations 5g(t) of g{t) vanishing at the endpoints. 



ii g(t) satisfies the Euler-Lagrange equations for L ao on G. 



iii The constrained variational principle 




(8.62) 



holds on q x V* , using variations of the form 



(8.63) 



where rj(t) £ q vanishes at the endpoints. 



iv The Euler-P oincare equations hold on g x V* 




(8.64) 
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Right Representation and Left Invariant Lagrangian. The set up is as fol- 
lows: 

• There is a right representation of Lie group G on the vector space V and G 
acts in the natural way on the left on TG x V*: h(v g , a) = (hv g , ah^ 1 ). 

• Assume that the function L : TG x V* —> R is left G-invariant. 

• In particular, if ao £ V*, define the Lagrangian L ao : TG — > R by L ao (v g ) = 
L(v g , do). Then L ao is left invariant under the lift to TG of the left action of 
G ao on G. 

• Left G-invariance of L permits us to define / : g x V* — > R by 



Conversely, this relation defines for any Z : q x V* — > R a left G-invariant 
function L : TG xP-tl. 

• For a curve G G, let £(t) := c/(i) — 1 _gr(t) and define the curve a(t) as 
the unique solution of the linear differential equation with time dependent 
coefficients a(t) = a(t)£(t) with initial condition a(0) = ao- The solution can 
be equivalently written as a(t) = a$g{t). 

Theorem 8.3 The following are equivalent: 
i Hamilton's variational principle 



Kg 1 v g ,a g) = L(v g ,a ). 




(8.65) 



holds, for variations 5g(t) of g{t) vanishing at the endpoints. 



ii 



g(t) satisfies the Euler-Lagrange equations for L ao on G. 



iii 



The constrained variational principle 




(8.66) 



holds on g x V* , using variations of the form 



6£ = ri + [£,rj\, 6a = ar], 



(8.67) 



where rj(t) £ Q vanishes at the endpoints. 



iv The Euler-Poincare equations hold on q x V* 



d 51 „ 51 51 



(8.68) 
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